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ABSTRACT 


This report covers the design of electric motors which fulfill the needs of Thrust Vector 
Control (TVQ actuators used in large rocket propelled launch vehicles. To accomplish this 
end the methodology of design is laid out in some detail. In addition a point design of a 
motor to fulfill the requirements of a certain actuator specified by MSFC is accomplished and 
reported upon. In the course of this design great stress has been placed on ridding the 
actuator of internally generated heat. To conduct the heat out of the motor use is made of the 
unique properties of the in house MSFC designed driving electronics. This property is that as 
long as they are operated in a quasi-linear mode the electronics nullify the effects of armature 
inductance as far as the phase of the armature current versus the rotor position is concerned. 
Actually the additional inductance due to the extended end turns in this design is of benefit 
because in the shorted armature failure mode the armature current in the fault (caused by the 
rotor flux sweeping past the armature) is diminished at a given rotor speed when compared to 
a more conventional motor with lower inductance. The magnetic circuit is analyzed using 
electromagnetic finite element methods. 
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STATEMENT OF RESEARCH PROBLEM AREA 


The specific task undertaken concerned itself with the design of electric motors for use with 
electromechanical actuators (EMA) used for thrust vector control. The outcome of the task 
was to provide insight into the design of motors in general, to provide insight into and 
solutions to problems unique to this EMA application and to provide a point design for an in- 
house MSFC EMA design. 

It was stipulated that brushless DC motors were the baseline technology. The motors were to 
be driven by MSFC designed electronics (a fortuitous choice as it turned out). 

The physical envelope into which the motors were to fit was determined by an existing 
MSFC actuator design. As it turned out this forced several features of the design. See 
discussion of results. 
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ANALYTICAL BACKGROUND 


All the phenomena of electromagnetic may be traced to Maxwell's equations. Until recently 
their direct use in low frequency analysis situations was limited because of the difficulty of 
solving them to obtain engineering results. This situation has changed somewhat with the 
advent of finite element approaches and the companion use of the computing power which 
has become available^). 


Maxwell's equations are listed below for convenience ( 2 ). 


VxH = i + 


ap 

at 


V x E = 


dB 

at 


Where £1 is the magnetic intensity, £ the electric vector, £ magnetic flux density, £ the 
electric flux density and i the conduction current The term 


dD 

at 


is the displacement current and at low frequencies such as are encountered in motor work 
may be assumed to be negligible (zero). 

In as much as the magnetic field is globally solenoidal (sourceless) the divergence of £ is 
zero and it is valuable as well as possible to find a magnetic potential vector, A- This is done 
by noting that £ is a vector and thus may be considered to be the curl of another vector 

(A)* 2 ). 

B = VXA 
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This development is particularly relevant to this work because the finite element program 
ANSYS produces A values as part of its operation. A may in turn be used to evaluate the 
back emf in the motor at any given speed. This development is as follows. From Faraday's 
Law 


EMF = 


ao 

9t 


where the flux is given by (for a given surface) 


O = j B • da 

and the voltage between two points is given by the line integral 

v u = Jkds 

Substitution yields 


fEds = -^ = -J^=da* 
Jl dt dt ' 


Applying Stokes' Theorem from vector analysis converts the line integral on the left to a 
surface integral. 


J(VxE)-da = -J^-da 

at 


From this it is seen that 


V x E = 


as 

dt 
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Making use of the magnetic potential formulation (B = VxA) then 


So that 


VxE 

“ ” 9t 


E = 


8A 

dt 


In this work interest is in voltage induced along one side of a coil of wire placed in the motor 
stator slots. This defines the directions of interest Letting A be directed along one side of 
the coil (down the slot) the electric vector £ may be calculated provided the speed of the 
motor vector is specified (i.e. to make a numerical approximation to At) and A A is found 
from A in the wire or coil side for different angles between the stator and the rotor magnets. 
Multiplying |E| by the length of the axially directed coil side yields the voltage of one coil 
side. Thus the voltage around the coil is given by 

V = 2\E\tn 

where n is the number of turns in the coil, t is the length of a coil side in a slot and the 2 
accounts for the coil having two sides in which voltage is being induced. This approach is 
different from the more familiar expression 

e = bfv (3) 


where v is the relative velocity of the flux density B and the length l. This is a simplified 
expression in that B and l are assumed orthogonal, the more correct expression is 


e = (B-v)/<” 
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This latter formulation is sufficiently removed from the first enumerated approach using A to 
be used as a "sanity check” on the results of the approach using A (which should be more 
accurate if fi is estimated from an ANSYS flux plot). 


Likewise there is a dual approach to calculating the torque developed by an electric motor. 
The more accurate one used by ANSYS is to evaluate the torque by means of the Maxwell 
stress tensor. ( 4 ) (5) The form of it useful here is 


This formulation is the one leading to the conceptual device of visualizing elastic strings 
(rubber bands) with one end attached to the stator and the other to the rotor so that as the 
stator associated mmf moves it tends to drag the rotor with it. Put a little more formally the 
magnetic field is attributed to posses stress analogous to stress in a structural member and in 
each case is associated a force on the surrounding environment S imil ar to the voltage 
calculation there is a better known way in which to calculate force on a conductor and hence 
torque on a rotor. It is to use the formula (3) 


T = rF = Bilr 


where T, F, B and i (or 1, in which i flows) are all orthogonal. In using this formula it is 
assumed that the current conductors which are really in the slots, are uniformly distributed 
around the stator side of the air gap. As Liwschtz-Garik has noted experiment confirms the 
legitimacy of this procedure^ 7 ). 
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Although the type of motor under consideration is sometimes called a brushless dc motor 
(because of the obvious lack of brushes and the fact that its power source can ultimately be 
traced to a dc bus) it is in reality more properly thought of as a synchronous motor. There are 
two major differences between this motor and the classical synchronous motor. One 
difference is that the electrical supply is controllable both as to current magnitude and phase 
(relative to the rotor-stator instantaneous geometry) and its temporal frequency (and hence 
the synchronous speed of rotation). The other is that the flux due to the rotor comes from 
permanent magnets rather than from electromagnets. Any course about electrical machines 
develops the equivalent circuit of the synchronous machine. References are numerous, see 
for instance numbers six and seven. This is done on a per phrase basis under the assumption 



This equivalent circuit is of one phase, the a phase, and hence the various a subscripts. Lms 
is the usual magnetizing inductance (with the reactive VARS drawn from the a source) here 
shown nonlinear, Ws is the synchronous speed (rad/sec), Wm is the mechanical speed of the 
rotor, Ta the torque assigned to the a phase, Rs and Ls the stator resistance and inductance, 
Nse and Nre the stator and rotor effective number of turns and la and Va the a phase amps 
and volts. The rotor is characterized by the field resistance (Rf) and inductance (Lef) and 
coupled to the stator through the ideal transformer as shown. 
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Reference six, through a series of transformations particularizes this general equivalent 
circuit to that corresponding to a cylindrical permanent magnet rotor design. It is insightful 
to point out some interesting steps in evolving permanent magnet equivalent circuits. This 
development parallels that in reference six to where the interested reader is directed for a 
more detailed account 

Consider the permanent magnet - soft iron circuit shown below in which the mmf required to 
established flux in the soft iron is assumed to be negligible. 



The permanent magnet has a cross section of area Am and length lm. The next figure shows 
the relationship between the magnetic flux 0m and the mmf 3n due to the current i in the N 
turns. 
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This graphical relationship between flux and mmf is obtained by rescaling a B-H curve 
no rmall y furnished by the material supplier to the dimensions of the specified magnetic 
circuit It should be noted that point a is where the demagnetization curve and the recoil line 
part company and that if the material is driven to the left of point a permanent weakening of 
the magnet will occur. It should also be noted that as the temperature of the magnet rises 
such a point will shift to the right Thus as the magnet temperature rises it becomes 
progressively easier to weaken the magnet It has been reported that such behavior may 
already have been observed in some designs. For this reason the designs done herein 
intentionally oversize the magnet (compared to using the maximum energy design points 
usually recommended). 


It will be appreciated in the model that a perfect transformer is being modeled in that there is 
no leakage accounted for. In order to obtain a linear model operation must be restrained only 
to the linear recoil line. The flux present in the core is found on the straight line which may 
be expressed as 

0m = 3m ^ +0rl 


where 


0rl 

3e 


_ 1 _ 

Re 


this may be manipulated to 

Ro0m = 3m + Ro0r/ 
or 

Ro0m = 3m + 3o 
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From this follows the equivalent circuit 



This in turn may be converted into an equivalent electric circuit as follows 

i = - 
1 N 

e = N— 
dt 

3 = Rd> 


Substituting 



di 

dt 
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Thus an equivalent circuit of the form 


Xe 



emerges. This may then be substituted for the rotor circuit of the original equivalent circuit 
to yield a per phase steady state equivalent circuit as 




The relevance and importance of the foregoing is that there is firm ground to think of a 
multiphase motor on a per phase basis provided it is being operated in a balanced fashion. 
This fact is used in sizing the motor. 
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MOTOR DESIGN CONSIDERATIONS 


Initially it was thought that a 5 inch diameter motor would be allowable. Subsequently it was 
decided that 3 inch diameter motors would be needed to fit a four electric motor actuator that 
MSFC has designed in house. From those exercises it seems that the somewhat larger 
diameter motor would be preferable because it allows keeping the speed at rated power and 
the length to diameter ratio smaller than is the case with the smaller diameter motor. 

The specification called for the actuator to develop 30 horsepower under nominal conditions 
but be able to develop 45 horsepower peak under off nominal conditions. Probably a desire 
to develop 30 horsepower at a stroking velocity of 5 inches/second and a stall force of 60,000 
pounds gave rise to these specifications (30 horsepower at 5 in/sec requires 40,000 pounds of 
force). Thus designing each motor for 7.5 horsepower no minal with a surge possibility of 11 
horsepower seems to be reasonable. The gear train connecting the motor to the load is two 
pass. One pass is a spur gear ratio in the range from 1 to 10. The other pass is roller screw 
gear with a range up to 1 inch per revolution possibility. Previous work ( 4 ) has shown the 
desirability of matching the reflected inertia of the load to that of the motor to the extent 
possible while still meeting the stall force requirement. 

When assigning a power rating to an electric motor close attention needs to be paid to the 
methodology followed. Electrtic motors are basically power limited by heating. Heating of 
the insulation of the wire used in the coils and, in a permanent magnet design, by magnet 
heating are principal questions. Overheating the wire insulation will drastically shorten its 
life even if not driven to instantaneous failure. The permanent magnets, as already noted, 
will, even below their Curie temperature, be susceptible to permanent irrevocable weakening 
because of the shift with temperature of the point where the demagnetization curve and the 
recoil line diverge. 
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Thus various ways of stating power ratings are encountered (sometimes for the same 
machine). Terms such as continuous rating, peak rating, time limited rating (all these with 
side conditions on the ambient environment, duty cycle, type of housing etc.) are found. 
They are codified for industrial use in the NEMA (national electrical manufacturers 
association) standards. The question here is how to rate the motors used in 
electromechanical actuators (EMA) used for large launch vehicle thrust vector control 
(TVC). 


Inspection of the flight records of Shuttle flights discloses very little activity of the actuators. 
The identifiable active times are during die initial roll maneuver to obtain orbital azim uth and 
during the maximum dynamic pressure encounter (max Q) when the air loads on this 
asymmetrical vehicle peak. But none of these loading periods last very long. Indeed the 
whole powered booster flight is about two minutes (without on average much activity). Thus 
short of an unspecified anomaly something other than flight stress needs to be used. 
Experience has shown many times that test and checkout of equipment stresses it more than 
flight. Operation in a boatail of a stacked vehicle at KSC in August during a functional 
checkout might be a high temperature regime. Perhaps others exist also. The current 
hydraulic actuator specifications (8) call for "operating performance...during and after 
exposure to..." an ambient of -40 to +275*F (-40*C to + 135*Q. The low temperature end of 
the scale would seem to impinge on electronic and mechanical things rather than 
electromagnetic concerns. The upper temperature range, when one considers that it begins to 
approach the upper limit of allowable magnet temperature is a problem. The magnets are 
innermost in the motor material and while heating of the rotor will be minimized by every 
means (e.g. using a laminated rotor) still heat generated there must cause a sensibly higher 
temperature than the ambient in order to make its way out through the air gap (radiation, 
perhaps convection) and the shaft (through the bearings). MSFC (Mr. Martin, EP-64) has 
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indicated upon questioning that the upper temperature might be modified as it is derived 
from a standard aerospace specification. To grasp the problem more quantitatively an 
estimate of the heat transfer possible from the body of the actuator was made as detailed in 
Appendix A. This development assumed that the motor heat is to be dissipated from the aft 
body or tailstock end of the actuator and that the actuator may be approximated by an 
equivalent right circular cylinder. For worst case conditions it was also assumed that the 
cylinder was horizontal. The natural convection heat transfer coefficient was obtained from 
reference nine. The results of the analysis are presented in Appendix A as the figure on page 
37. They show that heat dissipation during a stalled or otherwise loaded condition (i.e. high 
current) will be a problem. It may be ameliorated by the use of fins, perhaps a brake with the 
motor unenergized during motionless periods (adds complexity) or in extremum use of 
cooling air from a fan or other source. Once again the necessity of providing a low thermal 
resistance path from the origin of the heat to the actuator body is made plain. Appendix B 
does show for a representative geometry that if the thermal conductivity of the end turn to 
copper heat sink can be made approximately comparable to that of one inch of the wire itself 
than long periods of stalled operation could be thermally tolerated. 

In appendix C the core loss analysis is detailed. It shows that even under worst case 
conditions the temperature rise in the stator iron is small (delta temperature of 10-11 *F above 
stator outer temperature). Thus this factor is not a serious temperature rise concern. 
Comparing the wire temperature analysis to this one shows clearly that if a way could be 
found to couple the wire thermally to the stator iron an appreciable amount of heat could be 
extracted from the wire (recall that in appendix B it was assumed that all the heat exited the 
copper wire via the end windings). This might be brought about if the wires could be potted 
into the slots with something like the Stycast material. This latter might present fabrication 
problems and would have to be investigated. 
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The lamination design was driven by two factors. One was the fixed outer diameter (OD) of 
three inches. The other was the choice of the number of poles. The choice of the number of 
poles for a synchronous motor intended for conventional service is straight forward. For a 
three phrase p-pole, synchronous machine the relationship between mechanical shaft speed 
and the angular frequency of the supply is 

wm = — ws rad/sec. 

P 

Thus if the supply is 60HZ the only useful speeds i.e. those at which mechanical power can 
be produced^) are the discrete values 3600 rpm, 1800 rpm, 1200 rpm etc. Therefore to 
choose the number of poles the load speed characteristic is considered and a straight forward 
choice made. No such situation obtains here. In effect the electronics powering the motor 
can (and do) produce a wide range of frequencies. Thus some other criterion was needed. 
The initial choice made was to optimize the various possibilities according to the maximum 
torque they produce when full pitched windings are used. As seen below this led to an initial 
choice of a 4 pole design. To do this a series of 3 inch OD laminations were designed 
according to the procedure outlined below. They were then modeled in ANSYS and torque 
versus torque angle data gathered. Laminations incorporating 8, 6, 4 and 2 poles were 
designed and so evaluated. As is seen in the table below the torque production peaked for 
the 4 pole machine. However when later the zero current or reluctance torques were 
evaluated the choice changed to that of the 6 pole motor. 
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Nur 

Poles 

uber of 

Slots 

Toi 

Max. 

3" 

ique 

No Current 

< 

Toi 

Max. 

5" 

T}ue 

No Current 

2 

3x2x3 

24 

78 

27.3 



4 

3x4x2 

24 

108 

26.4136 



6 

3x6x2 

36 

96 

0.79968 

163.722 

-1.89312 

8 

3x8x2 

48 

77 

-4.68328 


— 


The units of torque are Newton-Meters per meter length of the rotor. 


















All the lamina tions were designed basically the same way. With the OD specified a priori 
the consensus formulas from reference fifteen (section 9-73) were used. But first it was 
necessary to decide how many slots would be specified for each number of poles. Also 
rationale pro or con fractional slotting was considered. Fractional slotting is often used when 
a lamination design exists and it is desired to use these laminations for a new design for 
which they were not designed. Other reasons to use fractional pitch windings include noise, 
wave form and locking torques specially on machines of eight or more poles( 15 X This is 
because it behaves as a winding with many slots per pole per phase which reduces the 
distribution factors of the harmonics. Clearly if the number of slots per pole per phase 
became less than two a fractional slot design would be called for because a single slot per 
pole per phase will in all likelihood lead to large harmonics^ 16 ). Induction or synchronous 
motors have essentially the same stator laminations and are designed in most cases with 
integral slots. Here there are less than eight poles and more than one slot per pole per phase. 
Therefore integral slots were chosen as appropriate and consistant with the general motor 
design literature. This choice renders the machine electrically easy to balance, makes 
analysis simpler because only one pole of the machine need be modeled (with appropriate 
boundary conditions) whereas a non integral machine must have the full 360 mechanical 
degrees modeled. However if testing shows the desirability of non integral slots and 
fractional pitch windings they could be designed. Two samples of commercial motors 
examined do indeed have non integral slot designs. The consensus formulas (dimensions in 
inches) are (given OD) 

inside diameter of rotor, D =-— 1 -- 

L175 + (1.03/P) 


where P is the number of poles. 


The diameter of the bottom of the slots is 


Di=1.175D +0.647 
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and the tooth width 


tw = 


L35D 

s, 


where Si is the number of slots. For partially closed slots the pole tip to pole tip distance e* 
is given by 

e 1 = 0.0143D +0.0643 


The value resulting from these formulas are likely dimensions for trial designs. However 
they have proven to be good guides. During the course of this design effort the outer radius 
of the motor (i.e. the OD) was increased as a trial. The torque analysis showed that even at 
high flux levels little torque gain resulted although of course the back iron flux density 
decreased. A value of one millim eter (about 4 mils) was chosen as the air gap between rotor 
and stator. This may be a little too small when the mechanical design is executed. A study 
was done which showed the air gap flux density to be a weak function of the air radial gap as 
long as that gap is s mall compared to the air gap between the alternating magnet edges. In 
the figure below the distance between the magnets was approximately 7 mm. 



Mr 0«p ttMtft (IMI 
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It is seen that another one or two millimeters would not be too bad although the gap should 
be kept as small as possible to achieve maximum air gap flux density. 

As previously discussed it was decided to have an integral number of slots per pole per 
phase. Thus for the six pole machine use of the pole /slot formula 

No. of slots = No. of phases X No. of poles X No. of slots per pole per phase 
or here yields the number of slots as 

36 = 3x6x2 

Applying the consensus formulas yields 


D = 1.7473 
Di = 2.7000 
tw = 0.0655 
el =0.0893 

It will be seen from the lamination figure that these values were adhered to by and large. 
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LAYOUT OF A 
6 POLE 3' OD 
MOTOR 
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This same procedure was followed for 2,4 and 8 pole designs. This was required to model 
them in ANSYS for the torque, flux and back EMF evaluations. 


As previously mentioned when discussing torque optimization all the work so far has been 
done with full pitch windings. However it may be desirable to chord or short pitch the stator 
windings to minimize harmonics. With two slots per pole per phase it would be possible to 
short pitch the winding by 30 electrical degrees (one pole is 180 electrical degrees, thus 180 
divided by 6 yields 30 electrical degrees per slot). A 150* pitch yields pitch, distribution and 
combined factors for various harmonics as (?) 


Harmonic 

EjKh.fl5.Q>Slfl 

Distribution Factor 

KpKd 

Fundamental 

0.966 

0.966 

0.933 

Third 

0.707 

0.707 

0.499 

Fifth 

0.259 

0.259 

0.067 

Seventh 

-0.259 

0.259 

0.067 

Ninth 

-0.707 

0.707 

0.499 

Eleventh 

-0.966 

0.966 

0.933 


From the table it is seen that the desired fundamental is diminis hed only 6.7% where as the 
other harmonics are diminished quite a bit Connecting the motor in wye or star will 
eliminate the third harmonic at the terminals and as a delta connection would allow 
circulating third harmonic currents, a wye connection is chosen. The de tails of the winding 
are given in appendix D. The flux plot see p. 21, of the 4 pole design is interesting. As 
previously pointed out only one quadrant of a 4 pole motor needs to be modeled 
(360*/4=90‘). The flux densities are shown in various colors ranging from dark blue to red. 
The peak air gap flux is quite high (above 2.9 Tesla) and perhaps just as 
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4 Pole 3" Motor 































important is oriented rather close to the optimum angle of 45* as derivable from the Maxwell 
stress expression based torque equation (4). Both of these factors explain the high peak 
torque production. At this point it is instructive to compare the flux plot of the 8 pole 
machine to that of the 4 pole machine. The plot is given on p. 23. Note that the overall 
levels of flux density are lower than for the 4 pole case. Also notice that there is only one 
area of 2.38 Tesla flux density per pole here whereas there were four areas of nearly 3 Teslas 
flux density per pole in the 4 pole case. Also notice that while much of the flux density 
orientation in the four pole case approximated the optimum 45 degree angle little of the flux 
(at a lower density level) in the eight pole case favored the 45 degree optimum. This 
spreading of the flux amongst a number of teeth may explain why the reluctance torque of 
the 8 pole motor is so low. The next figure, on p. 24, is a plot of torque (in Newton- 
Meter/Meter of motor length). From this it is seen that torque remains within approximately 
2 N-M/M of the maximum for commutation angles of 20* of peak. Thus it is seen that to 
miniTni7g torque ripple and retain a maximum torque capability at all shaft positions the use 
of a revolver based rotor - stator position sensor is desired rather than the 6 segment 3 Hall 
effect 60 electrical degree switching sometimes used. 
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Stress Plot of an 8 Pole 
















In [3 7]:= 

T4 » {{15,25.952),{16,27.347},{17,25.992},{18,26.906}, 
{19,24.727},{20,24.685},{21,24.417},{22,25.725}, 
{23,25.065},{24,25.183},{25,24.337},{26,23.6}, 
{27,23.091},{28,22.5},{29,21.984},{30,20.555}, 
{31,20.973},{33,20.02},{35,19.303},{37,18.303}, 
{39,16.548},{40,15.31}} 

ListPlot[T4,PlotJoined->True,GridLines->Automatic, 
AxesLabel-> { "Torque Angle Degrees", 

"Torque N-M/M"},PlotRange->{0,28} ] 


Torque N-M/M 
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Knowing that a full slot of 8 MA/M^ current density produces about 96 N-M/M of torque the 
motor sizing may be completed. For the ensuing calculations a fill factor of 50% will be 
assumed which means the peak torque actually available will be 96/2=48 N-M/M. The 
problem now is to establish the length of the active rotor and the speed at which rated power 
is to be developed. To keep the rotor length as short as possible the highest rotation speed 
thought practical i.e. 20,000 rpm was chosen. The rotor length calculation, assuming 48 N- 
M/M of torque is available to produce 7.5 horsepower is a follows. 


48 


N-M 

M 


x 0.7376 = 35.4——— 
m 


_ lb-ft ,.„_lb-ft 

7.5horsepower x 550-= 4125- 

sec- HP sec 


20,000 rpm = 2094 rad/sec 


2094— x 35.4f —-) x £ (m) = 4125f• l ^—^] 
sec V. m ) v sec ) 


e = 


4125 

(35.4)(2094) 


= 0.05565meter 


or 

i = (0..05565X39.37) = 2.191inches 


The rotor force stress is calculated as follows. The equivalent maximum tangential force is 


(F • r)(Speed) = power 


F = 


power 


4125 


(radiusXspeed) j L75 

WijJ 


K2094) 
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F=27.02 lbs 


The rotor surface area is 

so the force density is 


This is in line with usual good industrial practice^). See discussion in Appendix I. 

As was noted earlier a late switch was made in the choice of 6 versus 4 poles. Evaluating the 
A values requires many ANSYS runs and hence was not repeated. The method remains the 
same for any number of poles. To determine the number of turns per coil per pole ANSYS 
was queried to produce the magnetic potential ( 4 ) at the various coil locations so that the 
procedure described early in this report could be applied. The value of A at each side of each 
coil location or slot was found for a variety of rotor - stator angles and then averaged. These 
data were entered into a spread sheet program to find the coil (one turn) induced voltage. 
The spread sheet results are displayed in Appendix E. These results were calculated 
assuming a 10,000 rpm rotor speed 4 poles and a 7 inch long rotor. Because the calculations 
are linear they may be rescaled to the 20,000 rpm rotor speed as (taking 9V per coil as the 
spread sheet maximum) 


Area = (1.75*) (2.2) = 12.1in2 


8 = 1 = ^ = 12334 
A 12.1 in 2 


9x 


20,000 7.78 
10,000 X 7 


_ _ volts, . 

= 20——(max.) 
coil 


In this design where the basic dc bus is 270 volts (which may sag to 240 volts) the terminal 
voltage value is not critical as long as the back emf of the motor never rises so high as to 
threaten the linear operation of the current source amplifier feeding die motor. It is estimated 
that a voltage no higher than 200-210 volts will achieve this. Other than this consideration 
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there is much more freedom to select maximum back emf i.e. the number of turns in a coil 
than is the usual case. In this completely balanced electrical design there is freedom to series 
or parallel the coils in a given phase provided the current rating of any given coil is not 
exceeded (and since usually all coils are wound with the same wire gage this means having 
the same current in all the coils). If the terminal voltage is applied to the wye connected 
motor it is expected that the new phase voltage will be on the order of 200-210 divided by 
V3 (i.e. 115 to 120 volts). At 20 volts per coil only one or two turn coils in series connection 
(4x20=80 or 4x20x2=160) of all coils would produce an emf less than that desired. If the 4 
coils were connected 2 by 2 in series and then in parallel 3 turn coils (3x[20+20]) would 
yield 120 volts. The third possibility is to parallel all the coils with 6 turns each to yield 120 
volts (6X20). Keep in mind that each turn of the coil can be made from wire which in its 
turn is made from a number of wire strands in parallel. Once again consulting with the 
winder would help to make the choice by seeing which option will lead to the highest fill 
factor. It is essential that after a motor rotor and stator have been fabricated and assembled a 
search coil of a few turns be placed in the stator slots and the actual voltage be measured. 
With that knowledge the final winding specification can be developed. 

In the course of developing the exact winding to be used it should be kept in mind that the 
mechanical power delivered to the shaft is proportional to the back emf times the current 
Thus as the various emf options are considered from the standpoint of winding ease or 
complexity there are direct implications to the current drawn. That is if the back emf is 
relatively low then the current drawn has to be relatively high so that their required product 
i.e. the power delivered) is maintained. Of course the converse is true for a high back emf 
case. 


In a precursor study (reference four) a great deal of effort was spent in studying the 
interrelationships between the various system inertias. The more or less invariant inertia 
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component is that of the load (here the reflected rocket motor mass) and the actuator load 
bearing parts. The variable factors to consider when designing the system include the values 
of the spur gear ratio and the roller screw ratio between translation and rotation (which 
among other things has a gear train like inertia reflection property) the pinion and mating 
gear inertia and the motor rotor inertia. As is shown in the previous final report if the two 
passes are selected to produce matching of the referred inertias (referred to any common 
point) then two optimum objectives are achieved. These are first die m ax im u m bandwidth of 
the actuator (not the motor motion but the output motion) and second the maximum energy 
transfer to the load. Proper matching does not however generally achieve the other actuator 
requirement i.e. achieving specified actuator stall force. It was further shown that given a 
matched situation there is an optimum way in which to "detune" the system so as to achieve 
the required stall force and yet minimize the loss of system bandwidth and power transfer 
efficiency. Clearly any energy stored in the motor rotor is not energy getting to the actuator 
and load combination. Thus minimizing the motor inertia subject to power and thermal 
considerations is highly desirable. However this must be done in the context of an end to 
end optimized design to achieve its full potential. Attention is invited to Appendix I as well 
as the former work's final report for further thoughts on the inertia question. 
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Conclusions and Recommendations 


An extensive array of analytical tools has been assembled to design and analyze electric 
motors. These have been applied to the point design of a brushless dc electric motor for use 
with large thrust vector control (TVC) electromechanical actuators (EMA). This array 
represents the current state of the art in motor design tools. This is exemplified by the use of 
finite element methods (specifically the ANSYS electromagnetic option 4.4A) to evaluate 
numerically the motor designs torque characteristics (e.g. torque versus torque angle), the 
optimum number of poles, the back emf, and to visualize the flux paths of the proposed 
designs. Use has also been made of a digitally implemented thermal network analyzer that 
was obtained from COSMIC (and which runs on the PC) to investigate heat transfer in the 
proposed motor. The designs of almost all electric motors are thermally driven. This design 
should promote heat removal from the motor to the greatest extent possible. To do this some 
traditional design principals were contravened with the knowledge that the MSFC electronics 
will control the current even if the motor inductance is higher than normal. 

It will be noted in the table listing torque for various motors and circumstances that a 5 inch 
diameter motor is indicated. This is because initially such a motor was thought suitable for 
the actuator. However in addition to the usual tests a number of other tests were made with 
this model. Basically these tests held the slot current density constant as the motor's slots 
were elongated. The results of these preliminary tests showed that placing three separate 
windings in the slots with the currents in phase created no particular problem. If the motor 
were to be powered from a conventional low impedance or voltage source of electrical 
power it would however be very difficult if not impossible to insure adequately-in-phase 
currents. This is because as the slots are elongated to accommodate the additional windings 
the per winding inductance would be markedly different one from another and thus when 
they are all driven from the same low impedance source the various winding currents would 
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not be in phase (as is required to produce maximum or additive torque). Furthermore no 
passive tuning would be able to compensate for this because of the changing frequency 
applied as motor speed is changed. Fortunately the MSFC amplifier removes this current 
phasing problem/concem because it incorporates sufficient feedback so that its output 
current is essentially independent of the load impedance. This latter happy circumstance 
opens the possibility of incorporating three electrically isolated windings in a single motor 
stator with one rotor (on a single shaft) with the resulting implications on possible 
redundancy schemes. 

This motor design incorporates an innovative heat dissipation scheme. This scheme makes 
use of the MSFC amplifier indifference to winding inductance. This indifference allows an 
increase in the length of the end windings of the motor so as to couple the stator windings 
thermally with a copper heat sink placed at each end of the motor. The design then is such as 
to invite the thermal coupling of the motor to the actuator so as to make use of the surface 
area of the actuator for heat dissipation to the environment 

A design study has also been performed which looks into the use of a Titanium retaining 
sleeve proposed to be fitted around the rotor magnet structure. The purpose of this sleeve is 
to assure that the magnets remain tightly fixed to the rotor at high rotor speeds. An ANSYS 
study of this design disclosed that the Young's Modulus given in the literature for the 
magnets is extremely low. A check with the manufacturer showed an error of a thousand in 
the value. The ANSYS program was run with the correct value of Young's Modulus and 
showed that a feasible heating shrinking scenario is possible for mounting such a ring. An 
interesting result was that as the diameter of the motor increased the heat shrinking technique 
became less feasible. 
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This study has resulted in the electromagnetic design of a motor. There now needs to be a 
Tyu.rhanirai design made, materials ordered and the design resulting therefrom mechanized. 
This will serve to explore the ramifications of this outside diameter and thermally driven 
design. If discrepancies are found between the laboratory results and the analysis they will 
need to be rationalized. 

When building the motor the mechanical parts (including the rotor magnets) should be 
assembled with but a few turn test coil of 6 slot pitch in the stator. The rotor should then be 
spun at various rpm's to validate the back emf analysis. After this has been done and any 
modifications to the winding design made the motor would then be wound in the final stator 
winding configuration. Of course then tests should be performed to determine the developed 
torque at various electric current levels. 

At the conclusion of this build and test cycle there will be a motor design tailored to the TVC 
application with special emphasis on the heat transfer questions. 

In the future if the opportunity should arise to configure a new design it is recommended that 
serious consideration be given to relaxing the small (here 3 inch) outer diameter constraint. 
This constraint in the present design has caused the loss of the opportunity to optimize the 
aspect ratio of the motor. A shorter but larger diameter motor would be worthy of further 
serious study and would allow the optimization procedure to include additional parameters 
(i.e. simultaneous inertia matching and stall torque considerations as shown in previous 
work) than does the present case. 

Computer sketches of the motor design are enclosed immediately below. 
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TYPICAL 2-D MOTOR DRAWING 
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Appendix A: 


Heat Loss From Actuator Body 

As mentioned in the main body of the text the following calculations were made to gain 
knowledge of the heat transfer possibilities from the actuator body to the ambient air (at one 
atmosphere of pressure). Almost certainly test and check out of actuators will take place at 
one atmosphere. Flight, of course, will rapidly decrease the ambient pressure and heat 
energy will have to be stored in the thermal capacity of the motor and actuator for the 
relatively short (less than 10 minutes) flight time. 

As previously mentioned the Rockwell Specification document W calls for an operating 
ambient of 275*F (135*C). Magnequench 30 (the permanent magnet material modeled in this 
study) can withstand 180*C (356*F) to perhaps a maximum of 200*C (392*F). Thus there is 
something on the order of 80*F to 100*F (26*C to 38’C) between the maximum ambient 
(275*F) and the magnet temperature. Because most of the heat will be generated in the stator 
this does not leave much temperature difference to promote the outflow of thermal energy. 
Clearly if definitive data as to the maximum actual ambient were available (and showed less 
than 275*F) the situation would be eased. 

The actuator body to which the motors are mounted is described in reference 10. For the 
purpose of this work the tailstock of the actuator was assumed to be the active heat loss part 
because the bolted flange between it and the headstock end of the actuator is normally a 
relatively poor heat transmission path. 


34 



The heat loss surface was taken to be a right circular cylinder of Aluminum 12 inches in 
diameter and 14 inches long (based on an estimate of the convecting surface of the housing; 
in principle fins could be added to aid convection). 


Assuming the actuator is small compared to its enclosure and that the enclosure is at the 
game, temperature as the ambient the simplified heat transfer equations may be written as (for 
radiation and convection). ( 9 ) (H) 


qr = SejA, (Tj 4 - T 2 4 ) 


qc = hc(T, — T 2 ) 


where 

qi = rate of heat flow (BTU/HR) 

he = Convective heat transfer coefficient (BTU/HR-*R) 

Ti = actuator temperature (*R) 

T 2 = ambient temperature (°R) 

8 = Stefan - Boltzman constant (0.174X 10*^BTU/FT^-HR-°R^) 
e = emissivity of Aluminum 

Combining radiation and convection yields. 

qc + qr = {0.27(T, -T 2 ) 5/4 +0.174X10‘*X10' 1 [T 1 4 -T 2 4 ]}a, 

A, = 2IW = = 3.67FT 2 

T 2 - 275 + 459.69 = 734.69* R. (maximum ambient temperature) 
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The convective heat transfer term uses the relationship from page 177 of reference nine 



where 


Do = diameter of cylinder (FT) 

The emissivity of Aluminum was taken from Table A-23 of reference nine and is 
representative of a number of Al uminum finishes in the 75-212*F range. 

Using the factor 0.293 to convert BTU/HR to watts and plotting the heat transfer for several 
combinations of Ti and T2 (see figure on page 37) it is seen that the housing would have to 
get fairly warm with respect to the ambient to dispose of appreciable heat From the figure, 
if the housing temperature T1 were at 100'F and the ambient was 75*F the heat loss would be 
approximately 40 watts and so forth. It is not unusual for one motor of this type to dissipate 
at least this much power at stall. Here there are four motors; so clearly peak power at stall 
could be a major problem if the motors are called upon to sustain a large stalled force against 
a load or are cycled during testing for prolonged periods at high force loading. 
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DeltaQ 


In [4 7]:- 

aa-0.27 ((tl+459.69)-(t2+459.69)) A (5/4) 
bb*0.174 10 A -9 ((tl+459.69) A 4-(t2+459.69) A -4) 

Atm*3.67 

QcQrF*(aa+bb) Araa 

Plot3D[0.293 QcQrF,{tl,80,340),(t2,75,275),AxaaLabal->{"Tl" 
, T2, "Watts"), FacaGrids-X^J.1 J 
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Appendix B: 


Considerations On Individual Wire Heating 

No mattpr how many or how few turns are called for or how many conductors will be used in 
parallel ("in hand") for the windings, wire of about AWG25 will be used. This is because of 
the requirement to fit the wire into the slots and to wind it on a long thin stator. The data in 
reference twelve indicates current density in the wire as follows. 

Standard Machine (Presumably Air Cooled) 

2 KA/in 2 to 6 KA/in 2 

or 3x 10 6 A/m 2 to9.3xl0 6 A/m 2 

Exotic Machines (Liquid Cooled) 

12 KA/in 2 

or 18.6 xlO 6 A/m 2 

Pulse Loading 
20 Win 2 

or 31 x lO^ AAn 2 

The maximum current recommended in reference thirteen for 18 AWG (the smallest listed) 
copper wire according to the National Electric Code is 5 amperes. This yields a current 
density of 

3.9 Win 2 

or 6.1 x 10^ AAn 2 
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which falls in the standard machine range. There are some exotic wire insulations such as 
ceramic coatings in existence but it is not clear if they are suitable for use in a motor this size. 
At 8 x 10 6 A/m 2 (5.1613 KA/in 2 ) single 25 AWG wire will carry up to 

3.173 x 1(H x 5161.29 = 1.64 Amps. 

It is co mm on design practice to assume that all the heat energy generated in the stator copper 
windings should be transported to and dissipated from the end windings. This is generally 
considered conservative by up to 20% of the heat to be dissipated. The rationale is that the 
coupling of the wire to the stator iron is nebulous at best (and high resistance) and that 
conduction through the iron is harder than conduction through the copper. Our laboratory 
experiment disclosed fairly good iron copper coupling, surprisingly, but no thermal 
resistance values were obtained (See Appendix G). Dissipating the heat in the end turns 
commonly leads to having a fan mounted on the shaft of the motor to blow air over the end 
windings, down the air gap and across the other end winding. This design expedient is not 
available here. This is a servo motor application i.e. the motor does not spend much time 
turning and may spend time stalled but producing torque. The design pursued here 
substitutes large copper heat sinks at both ends of the motor for the usual fan. However, it is 
still necessary to remove the heat from these copper sinks and to get it into the sinks in the 
first place. 

As a matter of fact approaches reminiscent of this one are found in other servo motors 
although they are not carried as far as is proposed in this design. 
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For the sake of preliminary design suppose the wire in the slot is 7 inches long, AWG 25 and 
that it is carrying 1.64 Amps. In addition suppose it is connected to the two end copper heat 
sinks by a thermal resistance at each end of the wire (see below). 


olV A 7 ' »wg25 /{far; 

-« Cu Wire 1.64 A I 


Assume the Cu heat sinks’ temperatures may be varied for analysis purposes. Divide the wire 
into 7 segments and represent them as lumped conductances of 


thermal conductance = thermal conductivity x cross sectional area 

length 


= 0.061 


BTU/SEC 2.51xHT*Ft 2 
_ FfF X 144 


(1 / 12)Ft 


Cond = 1.2795 x 10~ 6 BTU/SEC per inch 

*F 


A transient analysis also requires each node to have a thermal capacity. It is calculated as 


Capacity = Volume x Density x Specific heat 


Here the 


Volume = 


(E) 

2.517xl0~ <> \ 

U2 J 

144 J 


1.456xlO“ 7 Ft 3 
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so 


Capacity = (L456 x 10" 7 )(556.42X0.092) = 7.4564x10“* — 


The power dissipated in each one inch section of a wire was calculated for the worst case as 
follows. The resistance of one inch of wire is, at 20*C, given from references 2 and 13 as 

R20=2.7 x 10"3 ohms 

This value was adjusted to 180*C by the relationship 03) 

Rt = + 0.00393(t - 20)] 

where C = percent purity of the copper (used 97%) 
t = temperature in degrees Celsius 

Rt = 4.5 x 10"3 ohms 


Thus the power per inch is given by 


P = RI 2 = (4.5xl0“* )(L65) 2 = 0.0123watts 


or 


P = 


(0.0123)(5.689x 10 2 ) - 

-- % - L = LI62 x 10 5 

60 


BTU 

SEC 
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The accompanying data, see pgs. 45-63, show a number of things. As they should the wire 
node temperatures final values vary linearly with the end node temperatures, the wire 
temperature is a strong function of the wire to heat sink conductance (it was varied from 10 
times better to 10 times worse than that of a one inch piece of wire). 

Qualitatively then it is clearly shown that keeping the heat sink temperatures down (this 
design incorporates a sink at each end of the motor; not just at one end) and minimizing the 
thermal resistance between the wire and the heat sink are vital. 

An estimate of the wire to heat sink conductivity was made by assuming the following 
geometry. 

1" long piece of Cu AWG 25 wire encapsulated 

A wo 2$ 

Ctpptr 
Wirt 


* . 018*1 


Qj 


i/3r 


k 


Assume a reasonable model is (of potting compound) 


Area = 


0.018"xl" 
144in 2 / ft 2 


= L 25x10^ ft 2 




2.6xl0” 3 ft 
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Conductance per inch of wire length is found from: 
(2850 KT Stykast) 


Conductivity = 19.3 


BTU~in .. 1 Ft_ L61 BTU 
Ft 2 HR*F X 12 in FtHRT 


which compares to Cu as 


210 


BTU 

HRFfF 


Thus the potting to Cu conductivity ratio inverse is 


210 

L61 


=130.6 


Conductance = Conductivity 


Area 

Length 


Cond.=> 


BTU Ft 8 ' BTU 
THRFt Ft ~ *FHR 


here for one inch of wire 


Conduc tan ce = 


(L61)(L25x10“*) 

(2.6X10") 


= 0.077 


BTU 

*FHR 


or 


0.077 


BTU 

*FHR 


1 HR 
X 3600 SEC 


= 2.15xl0' s 


BTU 

TSec 
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Compare to conductance of 1" of Cu wire (AWG25) from end to end which is 


1.2795 x 10" 6 .5ZLL 
Sec’F 

Ratio 

CondofCmpd 2.15 xlO" 5 

-— r — lO. o 

Ccu 1.2795 xlO -6 

If this can be achieved one would expect Re in the model to be somewhere between equal to 
Ccu and Ccu/10. The following data were run on the thermal analysis program for this case. 
From the results one sees that with the heat sinks held constant at 200’F, the wire started at 
100’F, and the conductance from the wire to the heat sink set equal to that of 1 inch of 
AWG25 (end to end conductance) the final hot spot temperature of 272*F (133*C) was 
reached in approximately 5 minutes. 133’C is well within the insulation capacity of good 
grade magnet wire. 

KLnr 


Thermal Network Used In Wire Analysis 
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752 DATA 2, 
754 DATA 1. 
756 DATA 3, 
->58 DATA 1. 
7&C DATA 4, 
762 DATA 1. 
764 DATA 5, 
766 DATA 1. 
768 DATA 6, 
770 DATA 1. 
772 DATA 7, 
774 DATA 1. 
776 DATA 8, 
778 DATA 1. 
780 Q(2) = 
782 Q(3) = 
784 Q(4) = 
786 Q(5) = 

788 Q(6) = 
790 Q(7) = 
792 Q(8) - 
794 T(1) = 
796 T(9) = 


7.46E-6,2, 
28E-6,1,1. 
7.46E-6,2, 
28E—6,2,1. 
7.46E—6,2, 
28E-6, 3, 
7.4 6E-6,2, 
28E-6,4,1 
7.46E-6,2 
28E-6,5,1 
7.46E-6,2 
28E-6,6,1 
7.46E-6,2 
28E-6,7,1 
1.162E-05 
1.165E—05 
1.165E-05 
1.165E-05 
1.165E-05 
1.165E-05 
1.165E-05 
200 
200 


0 

28E—6,3 
0 

28E—6,4 
0 

1.28E-6, 5 
0 

28E-6,6 
0 

28E-6,7 
0 

28E-6,8 
0 

28E-6,9 
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STEADY 

STATE 

SIGNIFICANT DIGITS= 4 

RELAX 

ITER'NS= 14 


NODE 

TEMP - F 

TEMP - R 

Q-INPUT Q-NET 

1 


200.000 

659.690 

0.000 

0.000000000000000 

2 


231.772 

691.462 

0.000 

0.000000051235634 

3 


254.519 

714.209 

0.000 

0.000000125966594 

4 


268.183 

727.873 

0.000 

0.000000111325612 

5 


272.753 

732.443 

0.000 

0.000000082773582 

6 


268.223 

727.913 

0.000 

0.000000051015277 

7 


254.587 

714.277 

0.000 

0.000000025734055 

8 


231.846 

691.536 

0.000 

0.000000009502783 

9 


200.000 

659.690 

0.000 

0.000000000000000 


Elapsed computer time with printing= 3.333333333330302D-02 minutes 
Date is 02-25-1993 Time is 09:40:18 
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TIME= 

0.0000 RCMIN= 

2.91406 MIN NODE # 

8 ITERATIONS= 0 


NODE 

TEMP - F 

TEMP - R 


Q-INPUT Q-NET 

1 


200.000 

659.690 

0.000 

0.000000000000000 

2 


100.000 

559.690 

0.000 

0.000139619995025 

3 


100.000 

559.690 

0.000 

0.000011649999578 

4 


100.000 

559.690 

0.000 

0.000011649999578 

5 


100.000 

559.690 

0.000 

0.000011649999578 

6 


100.000 

559.690 

0.000 

0.000011649999578 

7 


100.000 

559.690 

0.000 

0.000011649999578 

8 


100.000 

559.690 

0.000 

0.000139650001074 

9 


200.000 

659.690 

0.000 

0.000000000000000 

TIME= 

40.0000 RCMIN= 

2.91406 MIN 

NODE # 8 ITERATIONS= 0 


NODE 

TEMP - F 

TEMP - R 


Q-INPUT Q-NET 

1 


200.000 

659.690 

0.000 

0.000000000000000 

2 


205.622 

665.312 

0.000 

0.000001922702722 

3 


207.179 

666.869 

0.000 

0.000016218846213 

4 


204.963 

664.653 

0.000 

0.000004641791293 

5 


205.749 

665.439 

0.000 

0.000022936137611 

6 


204.969 

664.659 

0.000 

0.000004641791293 

7 


207.191 

666.881 

0.000 

0.000016219080862 

8 


205.639 

665.329 

0.000 

0.000001922702722 

9 


200.000 

659.690 

0.000 

0.000000000000000 

TIME= 

80.0000 RCMIN= 

2.91406 MIN 

NODE # 8 ITERATIONS= 0 


NODE 

TEMP - F 

TEMP - R 


Q-INPUT Q-NET 

1 


200.000 

659.690 

0.000 

0.000000000000000 

2 


223.042 

682.732 

0.000 

0.000001318639647 

3 


238.489 

698.179 

0.000 

0.000004236555014 

4 


247.018 

706.708 

0.000 

0.000003183478157 

5 


250.029 

709.719 

0.000 

0.000005991394573 

6 


247.024 

706.714 

0.000 

0.000003183478611 

7 


238.501 

698.191 

0.000 

0.000004236556379 

8 


223.059 

682.749 

0.000 

0.000001318640102 

9 


200.000 

659.690 

0.000 

0.000000000000000 

TIME= 

120.0000 RCMIN= 

2.91406 MIN 

NODE # 8 ITERATIONS= 0 


NODE 

TEMP - F 

TEMP - R 


Q-INPUT Q-NET 

1 


200.000 

659.690 

0.000 

0.000000000000000 

2 


228.875 

688.565 

0.000 

0.000000540234169 

3 


249.143 

708.833 

0.000 

0.000001254034032 

4 


261.101 

720.791 

0.000 

0.000001304240641 

5 


265.095 

724.785 

0.000 

0.000001773471922 

6 


261.106 

720.796 

0.000 

0.000001304240641 

7 


249.155 

708.845 

0.000 

0.000001254034032 

8 


228.892 

688.582 

0.000 

0.000000540234169 

9 


200.000 

659.690 

0.000 

0.000000000000000 

TIME= 

160.0000 RCMIN= 

2.91406 MIN 

NODE # 8 ITERATIONS= 0 


NODE 

TEMP - F 

TEMP - R 


Q-INPUT Q-NET 

1 


200.000 

659.690 

0.000 

0.000000000000000 

2 


230.837 

690.527 

0.000 

0.000000195862242 

3 


252.751 

712.441 

0.000 

0.000000398260966 

4 


265.838 

725.528 

0.000 

0.000000472853259 

5 


270.197 

729.887 

0.000 

0.000000563226081 

6 


265.844 

725.534 

0.000 

0.000000472853259 

7 


252.762 

712.452 

0.000 

0.000000398260966 

8 


230.855 

690.545 

0.000 

0.000000195862242 

9 


200.000 

659.690 

0.000 

0.000000000000000 
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TIME= 

200.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.498 

691.188 


0.000 

0.000000068028399 

3 

253.970 

713.660 


0.000 

0.000000130866653 

4 

267.434 

727.124 


0.000 

0.000000164235075 

5 

271.922 

731.612 


0.000 

0.000000185073404 

6 

267.440 

727.130 


0.000 

0.000000164235075 

7 

253.982 

713.672 


0.000 

0.000000130866653 

8 

231.516 

691.206 


0.000 

0.000000068028399 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

240.0000 RCMIN 3 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.721 

691.411 


0.000 

0.000000023233687 

3 

254.382 

714.072 


0.000 

0.000000043664503 

4 

267.973 

727.663 


0.000 

0.000000056091082 

5 

272.504 

732.194 


0.000 

0.000000061750931 

6 

267.979 

727.669 


0.000 

0.000000056091082 

7 

254.394 

714.084 


0.000 

0.000000043664503 

8 

231.739 

691.429 


0.000 

0.000000023233687 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

280.0000 RCMIN= 

2.91406 

MIN 

NODE i 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.797 

691.487 


0.000 

0.000000007880508 

3 

254.521 

714.211 


0.000 

0.000000014665628 

4 

268.155 

727.845 


0.000 

0.000000019025229 

5 

272.701 

732.391 


0.000 

0.000000020740330 

6 

268.160 

727.850 


0.000 

0.000000019025229 

7 

254.533 

714.223 


0.000 

0.000000014665628 

8 

231.814 

691.504 


0.000 

0.000000007880508 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

320.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.822 

691.512 


0.000 

0.000000002665296 

3 

254.568 

714.258 


0.000 

0.000000004939654 

4 

268.216 

727.906 


0.000 

0.000000006434593 

5 

272.767 

732.457 


0.000 

0.000000006985725 

6 

268.222 

727.912 


0.000 

0.000000006434593 

7 

254.580 

714.270 


0.000 

0.000000004939654 

8 

231.840 

691.530 


0.000 

0.000000002665296 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

360.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.831 

691.521 


0.000 

0.000000000900358 

3 

254.584 

714.274 


0.000 

0.000000001665751 

4 

268.237 

727.927 


0.000 

0.000000002173655 

5 

272.789 

732.479 


0.000 

0.000000002355728 

6 

268.242 

727.932 


O.ooo 

0.000000002173655 

7 

254.595 

714.285 


0.000 

0.000000001665751 

8 

231.848 

691.538 


0.000 

0.000000000900358 

9 

200.000 

659.690 

48 

0.000 

0.000000000000000 



TIME= 

400.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.833 

691.523 


0.000 

0.000000000303994 

3 

254.589 

714.279 


0.000 

0.000000000562008 

4 

268.244 

727.934 


0.000 

0.000000000733907 

5 

272.797 

732.487 


0.000 

0.000000000794799 

6 

268.249 

727.939 


0.000 

0.000000000733907 

7 

254.601 

714.291 


0.000 

0.000000000562008 

8 

231.851 

691.541 


0.000 

0.000000000303994 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

440.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.834 

691.524 


0.000 

0.000000000102618 

3 

254.591 

714.281 


0.000 

0.000000000189656 

4 

268.246 

727.936 


0.000 

0.000000000247742 

5 

272.799 

732.489 


0.000 

0.000000000268214 

6 

268.252 

727.942 


0.000 

0.000000000247742 

7 

254.603 

714.293 


0.000 

0.000000000189656 

8 

231.852 

691.542 


0.000 

0.000000000102618 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

480.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000034637 

3 

254.591 

714.281 


0.000 

0.000000000064007 

4 

268.247 

727.937 


0.000 

0.000000000083622 

5 

272.800 

732.490 


0.000 

0.000000000090520 

6 

268.253 

727.943 


0.000 

0.000000000083622 

7 

254.603 

714.293 


0.000 

0.000000000064007 

8 

231.852 

691.542 


0.000 

0.000000000034637 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

520.0000 RCMIN 3 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000011691 

3 

254.592 

714.282 


0.000 

0.000000000021603 

4 

268.247 

727.937 


0.000 

0.000000000028224 

5 

272.801 

732.491 


0.000 

0.000000000030551 

6 

268.253 

727.943 


0.000 

0.000000000028224 

7 

254.603 

714.293 


0.000 

0.000000000021603 

8 

231.852 

691.542 


0.000 

0.000000000011691 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

560.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000003946 

3 

254.592 

714.282 


0.000 

0.000000000007291 

4 

268.247 

727.937 


0.000 

0.000000000009526 

5 

272.801 

732.491 


0.000 

0.000000000010311 

6 

268.253 

727.943 


0.000 

0.000000000009526 

7 

254.603 

714.293 


0.000 

0.000000000007291 

8 

231.853 

691.543 


0.000 

0.000000000003946 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME 3 

600.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000001332 

3 

254.592 

714.282 


0.000 

0.000000000002461 

4 

268.247 

727.937 


0.000 

0.000000000003215 

5 

272.801 

732.491 


0.000 

0.000000000003480 

6 

268.253 

727.943 


0.000 

0.000000000003215 

7 

254.604 

714.294 


0.000 

0.000000000002461 

8 

231.853 

691.543 


0.000 

0.000000000001332 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

640.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000449 

3 

254.592 

714.282 


0.000 

0.000000000000831 

4 

268.247 

727.937 


0.000 

0.000000000001085 

5 

272.801 

732.491 


0.000 

0.000000000001175 

6 

268.253 

727.943 


0.000 

0.000000000001085 

7 

254.604 

714.294 


0.000 

0.000000000000831 

8 

231.853 

691.543 


0.000 

0.000000000000449 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

680.0000 RCMIN 3 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000152 

3 

254.592 

714.282 


0.000 

0.000000000000280 

4 

268.247 

727.937 


0.000 

0.000000000000366 

5 

272.801 

732.491 


0.000 

0.000000000000396 

6 

268.253 

727.943 


0.000 

0.000000000000366 

7 

254.604 

714.294 


0.000 

0.000000000000280 

8 

231.853 

691.543 


0.000 

0.000000000000152 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

720.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000051 

3 

254.592 

714.282 


0.000 

0.000000000000095 

4 

268.247 

727.937 


0.000 

0.000000000000124 

5 

272.801 

732.491 


0.000 

0.000000000000134 

6 

268.253 

727.943 


0.000 

0.000000000000124 

7 

254.604 

714.294 


0.000 

0.000000000000095 

8 

231.853 

691.543 


0.000 

0.000000000000051 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

760.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000017 

3 

254.592 

714.282 


0.000 

0.000000000000032 

4 

268.247 

727.937 


0.000 

0.000000000000042 

5 

272.801 

732.491 


0.000 

0.000000000000045 

6 

268.253 

727.943 


0.000 

0.000000000000042 

7 

254.604 

714.294 


0.000 

0.000000000000032 

8 

231.853 

691.543 


0.000 

0.000000000000017 

9 

200.000 

659.690 


0.000 

0.000000000000000 


50 



TIME= 

800.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000006 

3 

254.592 

714.282 


0.000 

0.000000000000011 

4 

268.247 

727.937 


0.000 

0.000000000000014 

5 

272.801 

732.491 


0.000 

0.000000000000015 

6 

268.253 

727.943 


0.000 

0.000000000000014 

7 

254.604 

714.294 


0.000 

0.000000000000011 

8 

231.853 

691.543 


0.000 

0.000000000000006 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

840.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000002 

3 

254.592 

714.282 


0.000 

0.000000000000004 

4 

268.247 

727.937 


0.000 

0.000000000000005 

5 

272.801 

732.491 


0.000 

0.000000000000005 

6 

268.253 

727.943 


0.000 

0.000000000000005 

7 

254.604 

714.294 


0.000 

0.000000000000004 

8 

231.853 

691.543 


0.000 

0.000000000000002 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

880.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000001 

3 

254.592 

714.282 


0.000 

0.000000000000001 

4 

268.247 

727.937 


0.000 

0.000000000000002 

5 

272.801 

732.491 


0.000 

0.000000000000002 

6 

268.253 

727.943 


0.000 

0.000000000000002 

7 

254.604 

714.294 


0.000 

0.000000000000001 

8 

231.853 

691.543 


0.000 

0.000000000000001 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

920.0000 RCMIN 3 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000001 

5 

272.801 

732.491 


0.000 

0.000000000000001 

6 

268.253 

727.943 


0.000 

0.000000000000001 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

960.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME=10 00.0000 RCMIN= 

2.91406 

MIN 

NODE # 3 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q-: 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

L040.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

1080.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

1120.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

1160.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME= 

1200.0000 RCMIN™ 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME™: 

1240.0000 RCMIN™ 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=1280.0000 RCMIN™ 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=1320.0000 RCMIN™ 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q-: 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=1360.0000 RCMIN™ 

2.91406 

MIN 

NODE # 8 


NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


TEMP - F 

TEMP - R 


200.000 

659.690 

0.000 

231.835 

691.525 

0.000 

254.592 

714.282 

0.000 

268.247 

727.937 

0.000 

272.801 

732.491 

0.000 

268.253 

727.943 

0.000 

254.604 

714.294 

0.000 

231.853 

691.543 

0.000 

200.000 

659.690 

0.000 


ITERATIONS™ 0 
PUT Q-NET 

0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
-0.000000000000000 
0.000000000000000 

ITERATIONS™ 0 
PUT Q-NET 

0.000000000000000 
0.000000000000000 
0.000000000000000 
-0.000000000000000 
0.000000000000000 
- 0.000000000000000 
0 . 000000000000000 
-0 . 000000000000000 
0.000000000000000 

ITERATIONS™ 0 
PUT Q-NET 

0.000000000000000 
- 0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

ITERATIONS™ 0 
PUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

ITERATIONS™ 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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TIME=1400. 

0000 RCMIN 3 

2.91406 

MIN 

NODE # 8 

NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=1440. 

0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


1 

2 

3 

4 

5 

6 

7 

8 
9 


TEMP 

200.000 

231.835 

254.592 

268.247 

272.801 

268.253 

254.604 

231.853 

200.000 


- F 


TIME=1480.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TIME 3 1520.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN 3 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TIME=1560.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TEMP 
659.690 
691.525 
714.282 
727.937 
732.491 
727.943 
714.294 
691.543 
659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 


- R 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


ITERATIONS 3 0 
JT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0 . 000000000000000 
0.000000000000000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0 . 000000000000000 
0.000000000000000 
0 . 000000000000000 
0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


MIN NODE # 8 ITERATIONS 3 0 

- R Q-INPUT Q-NET 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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TIME= 

1600.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME= 

1640.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

time=: 

L680.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q-: 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=1720.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


1 

2 

3 

4 

5 

6 

7 

8 
9 


TEMP 

200.000 

231.835 

254.592 

268.247 

272.801 

268.253 

254.604 

231.853 

200.000 


- F 


TIME=1760.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 


- R 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


ITERATIONS= 0 
JT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

ITERATIONS= 0 
JT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

ITERATIONS= 0 
JT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


ITERATIONS= 0 
Q-INPUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS= 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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TIME*: 

L800.0000 RCMIN* 

2.91406 

MIN 

NODE # 8 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME* 

1840.0000 RCMIN* 

2.91406 

MIN 

NODE # 8 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time*: 

1880.0000 RCMIN* 

2.91406 

MIN 

NODE # 8 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME* 

1920.0000 RCMIN* 

2.91406 

MIN 

NODE # 8 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time*: 

1960.0000 RCMIN* 

2.91406 

MIN 

NODE # 8 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME=2000.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TIME=2040.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TIME=2080.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TIME=2120.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TIME=2160.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 
659.690 
691.525 
714.282 
727.937 
732.491 
727.943 
714.294 
691.543 
659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 


MIN NODE i 
- R 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


f 8 ITERATIONS= 0 

Q-INPUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0 . 000000000000000 
0 . 000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS= 0 
Q-INPUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0 . 000000000000000 
0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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TIME 3 2000.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

2040.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME=2 080.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME=2120.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME=2160.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME=2200. 

0000 RCMIN 3 

2.91406 

MIN NODE # 8 

ITERATIONS 3 0 

NODE TEMP - F 

TEMP 

- R Q- 

INPUT Q-NET 

1 

200.000 

659.690 

0.000 

0.000000000000000 

2 

231.835 

691.525 

0.000 

-0.000000000000000 

3 

254.592 

714.282 

0.000 

0.000000000000000 

4 

268.247 

727.937 

0.000 

0.000000000000000 

5 

272.801 

732.491 

0.000 

0.000000000000000 

6 

268.253 

727.943 

0.000 

0.000000000000000 

7 

254.604 

714.294 

0.000 

0.000000000000000 

8 

231.853 

691.543 

0.000 

0.000000000000000 

9 

200.000 

659.690 

0.000 

0.000000000000000 

TIME=2240. 

0000 RCMIN 3 

2.91406 

MIN NODE # 8 

ITERATIONS 3 0 


NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


TEMP 

200.000 

231.835 

254.592 

268.247 

272.801 

268.253 

254.604 

231.853 

200.000 


- F 


TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 


- R Q-INPUT Q-NET 

0.000 0.000000000000000 
-0.000000000000000 
0.000000000000000 


TIME=2280.0000 


RCMIN= 2.91406 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

MIN NODE # 8 



NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=2320.0000 RCMIN 3 

2.91406 

MIN 

NODE # 8 


NODE TEMP 
1 200.000 

2 231.835 

3 254.592 

4 268.247 

5 272.801 

6 268.253 

7 254.604 

8 231.853 

9 200.000 


- F 


TIME=2360.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN 3 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 

2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 


- R 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

ITERATIONS= 0 
JT Q-NET 

0.000000000000000 

- 0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS 3 0 
Q-INPUT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
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TIME 5 *: 

2400.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

time=: 

2440.0000 RCMIN 5 * 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=2480.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q- 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 

TIME=2520.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 


NODE TEMP - F 

TEMP 

- R 

Q-: 

1 

200.000 

659.690 


0.000 

2 

231.835 

691.525 


0.000 

3 

254.592 

714.282 


0.000 

4 

268.247 

727.937 


0.000 

5 

272.801 

732.491 


0.000 

6 

268.253 

727.943 


0.000 

7 

254.604 

714.294 


0.000 

8 

231.853 

691.543 


0.000 

9 

200.000 

659.690 


0.000 


ITERATIONS= 0 
JT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

ITERATIONS* 5 0 
JT Q-NET 

0.000000000000000 
-0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

ITERATIONS 5 * 0 
JT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

ITERATIONS= 0 
JT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


TIME=2560.0000 
NODE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


RCMIN= 
TEMP - F 
200.000 
231.835 
254.592 
268.247 
272.801 
268.253 
254.604 
231.853 
200.000 


2.91406 

TEMP 

659.690 

691.525 

714.282 

727.937 

732.491 

727.943 

714.294 

691.543 

659.690 


MIN NODE # 8 
- R 

0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


ITERATIONS* 5 0 
Q-INPUT Q-NET 

0.000000000000000 

-0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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TIME=2600.0000 RCMIN= 

2.91406 

MIN 

NODE # 3 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

2640.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

2680.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=; 

2720.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME=2760.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME=: 

2800.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0 . 000000000000000 

2 

231.835 

691.525 


0.000 

- 0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0 . 000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

2840.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=; 

2880.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS^ 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

time=: 

2920.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 

TIME=2960.0000 RCMIN= 

2.91406 

MIN 

NODE # 8 

ITERATIONS= 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

200.000 

659.690 


0.000 

0.000000000000000 

2 

231.835 

691.525 


0.000 

-0.000000000000000 

3 

254.592 

714.282 


0.000 

0.000000000000000 

4 

268.247 

727.937 


0.000 

0.000000000000000 

5 

272.801 

732.491 


0.000 

0.000000000000000 

6 

268.253 

727.943 


0.000 

0.000000000000000 

7 

254.604 

714.294 


0.000 

0.000000000000000 

8 

231.853 

691.543 


0.000 

0.000000000000000 

9 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME=3000.0000 RCMIN= 

2.91406 

MIN NODE # 8 

ITERATIONS= 0 


NODE TEMP ■ 

- F 

TEMP 

- R Q-INPUT Q-NET 

1 

200.000 


659.690 

0.000 

0.000000000000000 

2 

231.835 


691.525 

0.000 

-0.000000000000000 

3 

254.592 


714.282 

0.000 

0.000000000000000 

4 

268.247 


727.937 

0.000 

0.000000000000000 

5 

272.801 


732.491 

0.000 

0.000000000000000 

6 

268.253 


727.943 

0.000 

0.000000000000000 

7 

254.604 


714.294 

0.000 

0.000000000000000 

8 

231.853 


691.543 

0.000 

0.000000000000000 

9 

200.000 


659.690 

1 0.000 

0.000000000000000 

Elapsed 

computer time 

with 

printing= 

1.383333333333326 

minutes 

Date is 

02-25-1993 

Time is 09:46: 

24 



752 DATA 2,7.46E-6,2,0 
754 DATA 1.28E-6,1,1.28E-6,3 
756 DATA 3,7.46E-6,2,0 
758 DATA 1.28E-6,2,1.28E-6,4 
760 DATA 4,7.46E-6,2,0 
762 DATA 1.28E-6, 3, 1.28E-6, 5 
764 DATA 5,7.46E-6,2,0 
766 DATA 1.28E-6,4,1.28E-6,6 
768 DATA 6,7.46E-6,2,0 
770 DATA 1.28E-6,5,1.28E-6,7 
772 DATA 7,7.46E-6,2,0 
774 DATA 1.28E-6,6,1.28E-6,8 
776 DATA 8,7.46E-6,2,0 
778 DATA 1.28E-6,7,1.28E-6,9 
780 Q(2) = 1.162E-05 
' Q(3) = 1.165E-05 
' - Q(4) = 1.165E-05 

786 Q(5) = 1.165E-05 
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Appendix C: 


Core or Iron Loss Considerations: 

Because this motor will be used in the servo application for TVC it will not in use experience 
much rotation and hence, in use, will not experience much core loss. However, during 
development, test and checkout it will be subjected to cyclic test, frequency response 
experiments and the like which will cause core loss. 

This analysis was carried out as follows. It was assumed that the motor was rotating at 
20,000 rpm and that this caused a sine wave of flux at a temporal frequency of 1333HZ. For 
the four pole motor, this should be a worst case by far. The separation of hysteresis and eddy 
current losses is difficult and not often attempted. Instead manufactures’ data are used. 
These data are derived in turn from test. Data from reference fourteen was cross checked 
with other data in reference six. From these references it was estimated that the core loss of 
this material at 1330HZ would be 45.35 watts per pound of material. The final four pole 
configuration was then divided into 24 identical segments. Each segment consisted of one 
tooth and the adjoining back iron half way to each of the adjacent teeth. A planar analysis 
was performed under the assumption that all heat flow of consequence will be radial because 
of the deliberately high axial thermal resistance in the stator stack (put there also to inhibit 
eddy current losses, which vary approximately as the thickness of the lamination squared). 
The use of the one twenty fourth segment is justified on the basis of symmetry i.e. all the 
segments are heating evenly and hence all heat flow must be radial. The segment was 
divided into six pieces and a lumped thermal network constructed as shown below. 
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2.272 2.615 2.97 4.444 4.18 

xlCT-6 xl 0*-6 xl 0*-6 xlCT-6 xio*-6 



Core Loss Thermal Network 


In the network above the initial dimensions are in millimeters but these were converted to 
feet in the calculations and the dimensions of the thermal conductance are 

BTU 

secT 

the dimensions of the thermal capacity are 

BTU 

•F 

and the entering power in 

BTU 

sec 

All temperatures are Fahrenheit From the diagrams it is seen that the outermost part of the 
stator was assumed to be at 200*F for the analysis. The rest of this appendix contains the 
network analyzer results and the input file used to characterize the network (see last sheet, 
lines 752 to 788). Both steady state and transient numerical analyses were performed. 
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STEADY 

STATE 

SIGNIFICANT DIGITS* 4 

RELAX 

ITER'NS* 12 


NODE 

temp 

- F TEMP - R 

Q-INPUT Q-NET 

1 


211.745 

671.435 

0.000 

-0.000000030378878 

2 


210.182 

669.872 

0.000 

-0.000000033029913 

3 


207.241 

666.931 

0.000 

-0.000000038733759 

4 


201.958 

661.648 

0.000 

0.000000007822710 

5 


202.473 

662.163 

0.000 

-0.000000014680492 

6 


201.958 

661.648 

0.000 

0.000000007822710 

7 


200.000 

659.690 

0.000 

0.000000000000000 


Elapsed computer time with printing* 3.333333333333144D-02 minutes 
Date is 02-26-1993 Time is 03:58:50 
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TIME= 

0.0000 RCMIN* 

0.97017 

MIN NODE # 

5 ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 


Q-INPUT Q-NET 

1 

100.000 

559.690 


0.000 

0.000002272000074 

2 

100.000 

559.690 


0.000 

0.000002615000085 

3 

100.000 

559.690 


0.000 

0.000002969999969 

4 

100.000 

559.690 


0.000 

0.000319479993777 

5 

100.000 

559.690 


0.000 

0.000340439990396 

6 

100.000 

559.690 


0.000 

0.000319479993777 

7 

200.000 

659.690 


0.000 

0.000000000000000 

TIME* 

40.0000 RCMIN= 

0.97017 

MIN 

NODE i 

5 5 ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 


Q-INPUT Q-NET 

1 

181.143 

640.833 


0-. 900 

0.000007742999514 

2 

184.844 

644.534 


0.000 

0.000007898842341 

3 

191.201 

650.891 


0.000 

0.000005324476206 

4 

199.929 

659.619 


0.000 

0.000000949694368 

5 

199.034 

658.724 


0.000 

0.000001700159146 

6 

199.929 

659.619 


0.000 

0.000000949694368 

7 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

80.0000 RCMIN= 

0.97017 

MIN 

NODE i 

5 5 ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 


Q-INPUT Q-NET 

1 

205.635 

665.325 


0.000 

0.000001540642529 

2 

205.125 

664.815 


0.000 

0.000001571254188 

3 

204.046 

663.736 


0.000 

0.000001058745056 

4 

201.554 

661.244 


0.000 

0.000000188528631 

5 

201.789 

661.479 


0.000 

0.000000337766863 

6 

201.554 

661.244 


0.000 

0.000000188528631 

7 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME= 

120.0000 RCMIN 3 

0.97017 

MIN 

NODE # 5 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

210.508 

670.198 


0.000 

0.000000306474334 

2 

209.160 

668.850 


0.000 

0.000000312563827 

3 

206.601 

666.291 


0.000 

0.000000210612242 

4 

201.878 

661.568 


0.000 

0.000000037503302 

5 

202.337 

662.027 


0.000 

0.000000067190712 

6 

201.878 

661.568 


0.000 

0.000000037503302 

7 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

160.0000 RCMIN= 

0.97017 

MIN 

NODE # 5 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q" 

INPUT Q-NET 

1 

211.477 

671.167 


o-.eoo 

0.000000060965810 

2 

209.962 

669.652 


0.000 

0.000000062177165 

3 

207.109 

666.799 


0.000 

0.000000041896318 

4 

201.942 

661.632 


0.000 

0.000000007460393 

5 

202.446 

662.136 


0.000 

0.000000013366001 

6 

201.942 

661.632 


0.000 

0.000000007460393 

7 

200.000 

659.690 


0.000 

0.000000000000000 

TIME= 

200.0000 RCMIN 3 

0.97017 

MIN 

NODE # 5 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

211.670 

671.360 


0.000 

0.000000012127704 

2 

210.122 

669.812 


0.000 

0.000000012368675 

3 

207.210 

666.900 


0.000 

0.000000008334280 

4 

201.955 

661.645 


0.000 

0.000000001484069 

5 

202.468 

662.158 


0.000 

0.000000002658849 

6 

201.955 

661.645 


0.000 

0.000000001484069 

7 

200.000 

659.690 


0.000 

0.000000000000000 

rlME 3 

240.0000 RCMIN= 

0.97017 

MIN 

NODE # 5 

ITERATIONS 3 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

211.708 

671.398 


0.000 

0.000000002412520 

2 

210.154 

669.844 


0.000 

0.000000002460455 

3 

207.230 

666.920 


0.000 

0.000000001657908 

4 

201.957 

661.647 


0.000 

0.000000000295220 

5 

202.472 

662.162 


0.000 

0.000000000528915 

6 

201.957 

661.647 


0.000 

0.000000000295220 

7 

200.000 

659.690 


0.000 

0.000000000000000 
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TIME* 

280.0000 RCMIN* 

0.97017 

MIN 

NODE # 5 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

211.716 

671.406 


0.000 

0.000000000479914 

2 

210.160 

669.850 


0.000 

0.000000000489449 

3 

207.234 

666.924 


0.000 

0.000000000329801 

4 

201.958 

661.648 


0.000 

0.000000000058727 

5 

202.473 

662.163 


0.000 

0.000000000105215 

6 

201.958 

661.648 


0.000 

0.000000000058727 

7 

200.000 

659.690 


0.000 

0.000000000000000 

TIME* 

320.0000 RCMIN* 

0.97017 

MIN 

NODE # 5 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

211.717 

671.407 


o.-ooo 

0.000000000095467 

2 

210.161 

669.851 


0.000 

0.000000000097364 

3 

207.235 

666.925 


0.000 

0.000000000065606 

4 

201.958 

661.648 


0.000 

0.000000000011682 

5 

202.473 

662.163 


0.000 

0.000000000020930 

6 

201.958 

661.648 


0.000 

0.000000000011682 

7 

200.000 

659.690 


0.000 

0.000000000000000 

TIME* 

360.0000 RCMIN* 

0.97017 

MIN 

NODE # 5 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

211.718 

671.408 


0.000 

0.000000000018991 

2 

210.162 

669.852 


0.000 

0.000000000019368 

3 

207.235 

666.925 


0.000 

0.000000000013051 

4 

201.958 

661.648 


0.000 

0.000000000002324 

5 

202.473 

662.163 


0.000 

0.000000000004164 

6 

201.958 

661.648 


0.000 

0.000000000002324 

7 

200.000 

659.690 


0.000 

0.000000000000000 

XIME* 

400.0000 RCMIN* 

0.97017 

MIN 

NODE # 5 

ITERATIONS* 0 


NODE TEMP - F 

TEMP 

- R 

Q- 

INPUT Q-NET 

1 

211.718 

671.408 


0.000 

0.000000000003778 

2 

210.162 

669.852 


0.000 

0.000000000003853 

3 

207.235 

666.925 


0.000 

0.000000000002596 

4 

201.958 

661.648 


0.000 

0.000000000000462 

5 

202.473 

662.163 


0.000 

0.000000000000828 

6 

201.958 

661.648 


0.000 

0.000000000000462 

7 

200.000 

659.690 


0.000 

0.000000000000000 


Elapsed computer time with printing* .8833333333333258 minutes 
Date is 02-26-1993 Time is 04:03:12 
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752 

DATA 

1,6.34IE—6,1,0 

754 

DATA 

1.4 6E-6,2 

756 

DATA 

2,7.81E—6,2,0 

758 

DATA 

1.46E-6,1,1.67E-6,3 

7C0 

DATA 

3, 8.31E-6,2,0 

762 

DATA 

1.67E-6,2,1.65E-6,5 

764 

DATA 

5,1.236E—5,4,0 

766 

DATA 

1.65E-6,3,3.865E-6,4,3.36E-6,7,3.865E-6,6 

768 

DATA 

4,1.17E—5,2,0 

770 

DATA 

3^8-65E--6, S^.aSSE-e^ 

772 

DATA 

6, 1.17E-5,2,0 

774 

DATA 

3.865E-6,5,3.153E-6,7 

776 

T (7) 

= 200 

778 

Qd) 

= 2.272E-06 

780 

Q(2) 

= 2.615E-06 

782 

Q(3) 

= 2.97E-06 

784 

Q(5) 

= 4.44E-06 

786 

Q(4) 

= 4.18E-06 

788 

Q(6) 

= 4.18E-06 
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Appendix D: 


Winding Considerations 

There are a n umb er of ways to depict winding information. The one preferred by this author 
is used below. It is assumed that a double layer i.e. two coils per slot winding is used. 
Because the lap and wave windings are electrically identical the choice between them is best 
left to the particular winding shop. 

First lay out a series of slots according to the number decided upon. In the bottom of the slot 
diagram fill the slots according to the slot star with the rotation CCW or in the positive angle 
direction. 



If a full pitch winding were being used both the top and the bottom of the slot would be 
carrying the same phase current If as is the case shown here on the next page the coils are to 
be short pitched by one slot they will span 5 slots rather than the 6 slots required for the full 
pitch case. This may easily be depicted by sliding the top row of conductors one slot to the 
left from its position in the full pitch case. The resulting diagram is shown on the next page. 
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The mmf resulting from such a winding may be depicted by assigning current magnitudes to 
the various phases (these magnitudes correspond to balanced 3 phase values e.g. A+=l, B+=- 
1/2, C+=-l/2, A-=-l, B-=l/2, C-=l/2). Then treat each slot as if the total current in that slot 
is the weight of an impulse and plot these impulses in the same position as each slot Next 
integrate spatially all the impulses and plot the results after having made the average of the 
mmf zero. A computer program was written to do this. It is written in C and produces 
Mathcad plot files. Thus a C compiler and Mathcad are needed to run it A listing is 
included below and a sample output is included. 
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/* This program calculates the MMF for an integer motor winding 

and creates files for the Mathcad program wind.mcd. 

Variable Definition: 

idisp = displacement of upper windings from lower 
winding, idisp 3 0 implies a full pitch 
winding, idisp 3 1 implies a winding short 
pitched by one slot, 
mmf 3 the magnetomotive force, 
iap * current in the the A+ phase 
icm * current in the C- phase 

ibp * current in the B+ phase 

iam 3 current in the A- phase 

icp = current in the C+ phase 

ibm * current in the B- phase 

phi * the phase angle of the A+ phase in degrees 
ampl = the amplitude of the currents 
current[i] 3 the magnitude of the ith current 
nslots 3 number of slots 
nspp 3 number of slots per pole 
nspppp 3 number of slots per pole per phase 
npoles 3 number of poles 
deg 3 conversion from degrees to radians 
ncases 3 number of cases; program allows calculation of 
up to 5 mmf's. 
icase 3 case counter 

1 include <stdio.h> 

# :lude <math.h> 

#j...elude <conio.h> 

Idefine PI 3.14159265 

Idefine deg PI/180, 

void main(void) 

< 

int nslots, nspp, nspppp, islot, ib, it, ipole, iphase; 
int npoles, i, icase, ncases, idisp, ip; 

float current[7], phi, ampl, mmf[201], iap, icm, ibp; 

float iam, icp, ibm, ibot[201], itop[201]; 

FILE *inp, *outl, *out2, *out3, *out4, *out5; 

inp 3 fopen("c:\\fort\\winding.in","r")? 
outl 3 fopen("c:\\fort\\mafl.out", W w"); 
out2 3 fopen("c:\\fort\\mmf2.out","w"); 
out3 3 fopen("c:\\fort\\mmf3.out", ,, w"); 
out4 3 fopen( w c:\\fort\\mmf4.out w ,"w''); 
out5 3 fopen( "c:\\fort\\mmf5.out", V); 


fscanf(inp," %d %d %d\n",incases, inslots, inpoles); 

for(icase 3 1; icase <» ncases; icase++) 

{ 

fscanf(inp,"%f %f %d\n", 

iphi, iampl, iidisp); 

nspp 3 nslots/npoles; /* # of slots per pole */ 

nspppp 3 nspp/3; /* # of slots per pole per phase */ 

iap 3 ampl*cos(phi*deg); 
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icm = ampl*cos((phi + 60.)*deg); 
ibp - ampl*cos((phi + 120.)*deg); 
iam = -iap; 
icp * -icn; 
ibm ■ -ibp; 

current(l] ■ iap; 
current!2] - icm; 
current!3] - ibp; 
current!4] - ian; 
current!5] * icp; 
current!6] » ibn; 

for (ipole ■ 1; ipole <- npoles; ipole +* 2) 

( 

for(iphase = 1; iphase <* 3; iphase++) 

( 

for(islot = l; islot <= nspppp; islot++) 

{ 

• ib = (ipole - l)*nspppp*3 + (iphase - l)*nspppp + islot; 
ibot(ib] = current[iphase]; 

) 

} 


for (ipole * 2; ipole <= npoles; ipole +=2) 

{ 

for(iphase * 4; iphase <» 6; iphase++) 

{ 

for(islot =1; islot <« nspppp; islot++) 

( 

■feib = (ipole - l)*nspppp*3 + (iphase - 4)*nspppp + islot; 
ibotfib] = currentIiphase]; 

> 

} 


for (ipole = 1; ipole <* npoles; ipole +**2) 

( 

for(iphase * 4; iphase <= 6; iphase++) 

( 

for(islot *1; islot <* nspppp; islot++) 

( 

ib - (ipole - l)*nspppp*3 + (iphase - 4)*nspppp + islot; 
ib - ib + nspp - idisp; 
if (ib > nslots) 

ib ■ ib - nslots; 
itopfib] * current[iphase]; 

} 

> 


for (ipole - 2; ipole <= npoles; ipole += 2) 

( 

for (iphase = l; iphase <= 3; iphase++) 

( 

for (islot =* 1; islot <* nspppp; islot++) 

{ 

ib - (ipole - l)*nspppp*3 + (iphase - l)*nspppp + islot; 
ib * ib + nspp - idisp; 

l' 
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} 


if (ib > nslots) 

ib ■ ib - nslots; 
itopfib] = current[iphase]; 

} 

) 

) 

mof CO] * 0.; 

for (islot * 1; islot <» nslots; islot++) 

{ 

mmf[islot] = mmf[islot - 1] + ibot[islot] + itop[islot]; 


switch (icase) 

( 

case 1: 

for(i = 0; i <=* nslots; i++) 

fprintf(outl," %f\n", mmf[i]); 
break; 
case 2: 

for(i =0; i <= nslots; i++) 

fprintf(out2," %f\n", mnf[i]); 
break; 
case 3: 

for(i = 0; i <* nslots; i*M>) 

fprintf(out3," %f\n", mmf[i])? 
break; 
case 4: 

for(i * 0; i <* nslots; i++) 

fprintf(out4," %f\n", mrf[i]); 
break; 
case 5: 

for(i =0; i <= nslots; i++) 

fprintf(out5," %f\n" # amf[i]); 
break; 

) 

) 

fcloseall(); 


5 24 4 

0 . 1 . 0 
30 . 1 . 0 
60 . 1 . 0 
90 . 1 . 0 

120. 1. 0 0.000000 
2.000000 
4.000000 
5.000000 
6.000000 
5 000000 
>00000 
2.000000 
0.000000 
>1.000000 
- 2.000000 



-1.000000 

0.000000 

2.000000 

4-000000 

>00000 

6.000000 

5.000000 

4.000000 

2.000000 

0.000000 

- 1.000000 

- 2.000000 

- 1.000000 

0.000000 




MMF = READPRN ( MMF5 OUT) 






8 



APPENDIX E 

Magnetic Potential Calculations 

This appendix contains the spread sheet which takes the averaged data from each side 
of the slots and for a 7 inch long rotor turning at 10,000 rpm produces the voltage that 
would be induced in a one turn full pitch coil of the 4 pole configuration. The 
maximum value of this voltage is used to determine the number of turns in the windings 
as explained in the main body of the text. 


81 



average (We/m Degree) Delta A (We/m De gree) [Electric Force (We/m si I Back EMF (Volts) 

_ 7-56E-03 _ -1.50E-05 _ -0.90423486 -0.321545916 

_ 3-49E-03 _ -5.05E-05 _ -3.0383145 -1.080424636 

_ -1-34E-03 _ 4.05E-08 _ 000243486 0.000865836 

_-5.82E03_ -7.85E-06 -0.47176164 -0.167758439 
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APPENDIX F 


Explanation of ANSYS44A Code for File MTRM15A (Six Pole Motor) 

The following program listing represents one sixth of the motor that was designed. 
This is a valid representation not only because the motor is symmetrical, but also because 
the analysis at all times spans one full pole. The electromagnetic flux that follows the air 
gap path between the stator and the magnet produces a torque, which effectively can be 
calculated by performing a line integral over that path. By definition, magnetic flux 
follows a closed path, so on either side of the motor, the magnetic flux should be equal 
and of opposite sign, therefore negating any contributions of flux on the sides of a pole to 
torque. The torque then, that a motor experiences may be calculated by summing up the 
line integral over a path in the air gap for a single pole to get a torque per pole value, and 
for a total torque value, multiplying the unit torque by the number of poles. 

ANSYS44A was instrumental in the simulation of this motor because of its ability 
to handle the nonlinear material properties inherent in magnetic materials. The code 
breaks up the design in several phases. The first of these phases is to define the material 
types by assigning numbers to materials and properties to the material numbers to be used 
in the analysis. Material one is defined to be air, and material two is the rotor iron 
(Carpenter HI-MU 49). When considering the rotor iron and also the magnets (materials 
3 and 4, which represent magnetization outward and inward respectively, 
MAGNEQUENCH MQ3-30H), the nonlinear properties of the materials must be 
included. This is accomplished by providing a BH curve for ANSYS44A to use, for each 
material posessing magnetic properties. In addition to these properties, for the magnets 
coercivity must be entered. For all of the material whether linear or nonlinear, the density 
of each material must be entered. A different material number was assigned to each phase 
of the windings: A+ A-, B+ B-, and C+ C-. 
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The next phase in this model was to create the nodes for the motor (see 
accompanying figure). This is accomplished by setting nodal limits. The radii from the 
center of the shaft to a material or shape change determine the limits. RO is the radius of 
the rotor shaft. R1 is the radius to the bottom of the magnets. R2 is the radius to the 
bottom of the air gap. R3 is the radius to the bottom of the shoes. R4 is the radius to the 
top of the shoes. R5 is the radius to the top of the slots. Finally, R6 is the radius to the 
outer motor. Nodes between the radii limits are generated from RO to Rl, R1 to R2, etc., 
such that there are a total of 31 nodes from RO to R6. It is important to have several 
elements in the air gap in the radial direction so that torque can be calculated accurately. 
For this model, five nodes in the air gap providing 4 elements was deemed sufficient. 

Once the basic node pattern has been established along the radial direction, this 
pattern is copied along the theta direction for 60 degrees. Sixty degrees is calculated by 
dividing the number of elements per magnet, plus the air gap elements between magnets 
by the number of poles. The width of each element is one degree and the motor spans 360 
degrees. Therefore, because this is a six pole motor simulation, 60 degrees is the sector of 
the motor over which the calculation is being done. The number of teeth per sector is 
found by dividing the number of elements per sector by the number of degrees per tooth 
(60/15), four slots per sector. The element type used for this model is the quadrilateral. 
Selecting a material number, for example, material two, is necessary to put the proper 
materials in their respective places, otherwise the default material number is one, which is 

air. Elements are formed by joining four nodes in this order 

• • 

3 2 

• • 

4 

1 

This pattern is repeated along the radial node path from element 1,2,33,32,1 to element 
30,31,62,61,30, by the command egen. Next, this swath of elements is copied along the 
theta direction for 60 degrees to generate all of the elements for the simulation All 
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elements were generated initially as rotor iron simply because the majority of the motor is 
composed of HOE-MU 49. The other parts of the motor such as: magnets, shoes, slots, and 
air gaps are formed by literally punching out a section of the HI-MU 49 and replacing it 
with those materials over the range of elements for which they span. 

The elements for the air gap were modified next by selecting the material number 
one, the nodal path that bounds the air gap in both the r and theta direction, and finally 
selecting all the elements within those bounds and modifying them to be composed of air. 
This procedure is also repeated for the magnets, shoes, and slots. The materials were 
assigned colors, to see if each material was in its correct place. 

For a finite element program solving the magnetostatic equations, the natural 
boundary condition (BC) is the zero Neumann boundary condition. This means that if a 
BC is not specified, the one applied is the zero normal derivative of the magnetic vector 
potential. This boundary condition is clearly not desired at the outer edge of the back iron 
because it implies that the flux is perpendicular to the outer edge of the stator. The 
appropriate BC is to fix the magnetic vector potential. Usually, the potential is set to 
zero. 

The appropriate BC at the inner shaft is not so easy to understand, but the 
condition that is commonly used is also the zero vector potential. Actually, the solution is 
not sensitive to the inner BC. 

The final phase was the actual analysis of the model. ANSYS44A uses the 
Maxwell Stress Tensor to determine the maximum torque in these electromagnetic 
models. The result of which shows that the maximum torque occurs in the air gap when 
the flux is at a maximum and is directed at 45 degrees. ANSYS44A calculates and stores 
the magnetic flux by first having been given a convergence criterion and a maximum 
number of load steps in which to reach this criterion by means of iteration. Once the 
iteration at which the solution converges has been isolated, that load step is written as the 
final load step and is saved to FILE12 for further analysis in the ANSYS44A post 
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processor (POST1). Once the solution is completed in the post processor, it can be used 
to calculate the torque by performing a line integral between the magnetic flux and the 
length vector over a specified nodal path in the air gap. 

These operations, once the initial model has been created can either be run 
interactively or in batch mode. In the interactive mode, a single simulation, it takes 
anywhere from 25 to 45 minutes. This is understandable because the magnetic flux must 
be evaluated at every node, and there are 1891 nodes and 1800 elements. In the batch 
mode, all of the steps are preprogrammed and all that is necessary to analyze the model is 
to submit the batch job. The IBM 3084 system takes about the same amount of time to 
run the batch job as the interactive, however, the batch job does not require constant 
human monitoring. Another advantage to the batch mode is that as many as four jobs at a 
time can be submitted. This is important because to show a fully developed torque curve, 
the model must be simulated at a range of rotor angles. Usually this requires at least 25 
angular displacements. 
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PROGRAM LISTING OF A 6 POLE 3" OD MOTOR 
DESIGN FOR ANSYS44A 


/CORE,5800000 

/CORANSYS REVISION 4.4A UP419 97 13.9922 11/14/1991 

1. /TITLE, EMA MOTOR ANALYSIS 

2. /PREP7 

3. KAN.-l 

4. KAY, 9,1 

3. ET,1,13^1 

6. KNL.1 

7. NLSIZE^OO 

8. /COM, 

9. /COM, MATERIAL 1 IS AIR 

10. /COM, 

1L NL>1,57,1 

12. /COM, 

13. /COM, MATERIAL 2 IS THE ROTOR IRON - CARPENTER HI-MU 49 

14. /COM, 

15. NLTABA5 

16. MP,DENS,2,8249 

17. NLX„1,10.,15.^20^25^30^35. 

18. NLX^.40,45,50,55,60,65 

19. NLX m 70,75,80,85^0^5 

20. NLX^lOO,795.8 

21. NLY^^.947,1.052,1.083,1.102,1.117,1.130 

22. NLY^l.140,1.147,1.155,1.164,1.172,1.181 

23. NLY^ 1.189,1.195,1.201,1.206,1.212,1.218 

24. NLY^l.222,1.5 

25. /SHOW,4208 

26. /COM, 

27. /COM, DEFINE MAGNET CURVE AS MATERIALS 3 & 4 

28. /COM, MAGNETIZATION OUTWARD FOR 3, INWARD FOR 4 

29. /COM, MATERIAL -MAGNEQUENCH MQ3-F30H 

30. /COM, 

31. CSYS,1 

32. NLTAB^,5 

33. MPDENS^.7550 

34. NIXU1800,827600 

35. NLY W 1^.06,1.13 

36. HC- -827600 

37. NL3,109JIC 

38. NL3,163,0 

39. NLTAB,4,5 

40. MP,DENS,4,7550 

41. NLX W 131800,827600 
4Z NLY W 1 W 0.06,1.13 

43. HC-827600 

44. NL*4,109,HC 

45. NL,4,163,0. 

46. /COM, 

47. /COM, MATERIALS 5 THROUGH 10 ARE A+, A-, B+, B-, 

48. /COM, C+, AND C- SOURCES RESPECTIVELY 

49. /COM, 

50. NL ,5,55,1 

51. NL ,6,55,1 

52. NL ,7,55,1 

53. NL ,8,55,1 

54. NL ,9,55,1 

55. NL ,10,55,1 

56. /COM, 

57. /COM, THET - THETA ORIENTATION OF ROTORS MEASURED TO LEADING 

58. /COM, EDGE OF OUTWARDLY MAGNETIZED MAGNET 

59. /COM, NTEE-NUMBER OF TEETH 

60. /COM, NMAG-NUMBER OF MAGNETS 

61. /COM, NTHE-NUMBER OF ELEMENTS IN THE THETA DIRECTION 

62. /COM, 

63. /COM, R0 - RADIUS IN METERS FROM CENTER TO OUTSIDE OF THE SHAFT 

64. /COM, Rl-RADIUS TO BASE OF MAGNET 

65. /COM, R2-RADIUS TO THE TOP OF THE MAGNET 

66. /COM, R3-RADIUS TO TOOTH SHOE AT AIR GAP 
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67. /COM, R4 - RADIUS TO END OF TOOTH AT THE SHOE 

68. /COM, R5 - RADIUS TO BACK IRON 

69. /COM, R6 - OUTER RADIUS OF THE MOTOR 

70. /COM, 

71. THET- -15. 

72. NTEE-36 

73. NMAG-6 

74. NTHE-360 

75. R0-0.005 

76. Rl-0.01118 

77. R2-0.02118 

78. R3 - 0.02218 

79. R4 - 0.02252 

80. R5 - 0.03429 

81. R6- 0.03810 

82. NR01-5 

83. NR12-9 

84. NR23-5 

85. NR34-2 

86. NR45-5 

87. NR56-4 

88. NTR1 - (NR01 + 1) 

89. NTR2 * NTR1 + NR12 

90. NTR3-NTR2 + NR23 

91. NTR4-NTR3 + NR34 

92. NTR5-NTR4 + NR45 

93. NTR6-NTR5 + NR56 

94. THMG-50 

95. /COM, 

96. /COM, THMA - NUMBER OF ELEMENTS PER MAGNET + AIR GAP BETWEEN 

97. /COM, NTPC - NUMBER OF TEETH PER SECTOR OVER WHICH THE CALCULATION 

98. /COM, IS DONE 

99. /COM, 

100. THMA-NTHE/NMAG 

101. THPT * NTHE/NTEE 

102. NTPS - THMA/THPT 

103. N,1,R0 

104. N,NTR1,R1 

105. FILL,1,NTR1 

106. N^moR2 

107. FILI^NTR1^TR2 

108. N,NTR3,R3 

109. FILUNTR2^rTR3 

110. N,NTR4,R4 

111. FTIJ^NTR3 y NTR4 

112. N,NTR5,R3 

113. FILL^TR4,NTR5 

114. N,NTR6,R6 

115. FILJUNTR5J4TR6 

116. DTHE - (360/NTHE) 

117. /COM, THMX - MAXIMUM THETA OVER WHICH WE ARE DOING THE CALCULATION 

118. THMX - NTHE/NMAG 

119. NELT - (THMX/DTHE) 

120. NGEN,(NELT+l)^rTR6,lXrR6,U043THE 

121. /COM, 

122. /COM, GENERATE ALL ELEMENTS AS ROTOR IRON 

123. /COM, 

124. NALLSEALL 

125. MAT.2 

126. N1 -1 

127. N2-2 

128. N3 - (NTR6 + 2) 

129. N4 - (NTR6 + 1) 

130. EJ41,N2^J3,N4 

131. ETR6 * (NTR6 -1) 

132. EGEN,ETR6,1,1,1 

133. EGEN,NELT,NTR6,1,ETR£,1 

134. /COM, 

135. /COM, MODIFY ELEMENTS FOR AIR GAP 

136. /COM, 

137. EPS-.00001 


92 



138. MAT.l 

139. NSEUX(R2-EPS)^3 

140. ENODE.1 

141. EMODIF,ALL 

142. EALLSNALL 

143. /COM, 

144. /COM, NOW DO THE MAGNETS - AIR 

145. /COM, 

146. MAT3 

147. NSELX(R1-EPS%R2 

148. ENODE.1 

149. EMODIFrALL 

150. MAT,1 

151. T1 - THMG + THET - EPS 

152. T2 -THMA + THET + EPS 

153. NRSEL*Y,T1,T2 

154. ENODE.1 SEMODIF, ALL 

155. NALL SEALL 

156. THEA-THMA-THMG 

157. Tl-THET-THEA-EPS 

158. T2 -THET + EPS 

159. NSELX,(R1-EPSXR2 

160. NRSEL.Y.T1.T2 

161. ENODE.1 SEMODIF^LL 

162. NALL SEALL 

163. NSELX(R1 -EPS),R2 

164. Tl - THET-THMA-THEA-EPS 

165. T2 - THET - THMA + EPS 

166. NRSEL.Y.T1.T2 

167. ENODE.1 SEMODIF,ALL 

168. NALL SEALL 

169. MAT y 4 

170. Tl - THET - THMA - EPS 

171. T2 - THET - THEA + EPS 

172. NSELXCRl - EPS),R2 SNRSEL,Y,T1,T2 

173. ENODE.1 SEMODEF.ALL 

174. NALL SEALL 

175. T1-THMA+THET-EPS 

176. T2 * THMA + THMG + THET + EPS 

177. NSELX,(R1-EPS%R2 SNRSEL,Y,T1,T2 

178. ENODE.1 SEMODIF^ALL 

179. NALL SEALL 

180. /COM, 

181. /COM, SHOES -AIR 

182. /COM, 

183. MAT.l 

184. /COM, THSH - 8.202 

185. /COM, THSA - 11.122 

186. THSH-8STHSA- 10 

187. THAR-THSA-THSH 

188. Tl-THSH/2 ST1-T1-EPS 

189. T2-T1 + THAR ST2-T2 +EPS ST2-T2 +EPS 

190. NSELX(R3-EPSXR4 SNRSEL*Y,T1,T2 SENODE.1 

191. EMODIF.ALL 

192. EALL SNALL 

193. Tl - T2 + THSH ST1 - Tl - EPS ST1 - Tl - EPS 

194. T2 - Tl + THAR ST2 - T2 + EPS ST2 - T2 + EPS 

195. NSELX(R3-EPSXR4 SNRSEL,Y,T1,T2 SENODE.1 SEMODIF^LL 

196. NALL SEALL 

197. Tl - T2 + THSH ST2 - Tl + THAR ST1 - Tl - EPS ST1 - Tl - EPS 

198. T2-T1 + THAR ST2-T2 +EPS ST2-T2 +EPS 

199. NSELX^R3-EPSXR4 SNRSEL,Y,T1,T2 SENODE.1 SEMODIF r ALL 

200. EALLSNALL 

201. Tl - T2 + THSH $T1 - Tl - EPS ST1 - Tl - EPS 

202. T2 - Tl + THAR JT2 « T2 + EPS ST2 - T2 + EPS 

203. NSEL,X(R3-EPSXR4 SNRSEL,Y,T1,T2 SENODE.1 SEMODIFRALL 

204. NALL SEALL 

205. Tl - T2 + THSH $T1 - Tl - EPS ST1 - Tl - EPS 

206. T2 - Tl + THAR ST2-T2 +EPS ST2-T2 +EPS 

207. NSELX(R3-EPS)J14 SNRSEL,Y,T1,T2 SENODE.1 SEMODIFRALL 

208. NALL SEALL 
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209. T1 - T2 + THSH $T1 - T1 - EPS ST1 « T1 - EPS 

210. T2 - T1 + THAR ST2 - T2 + EPS $T2 - T2 + EPS 

211. NSEUX(R3-EPS\R4 SNRSEL, Y,T 1,T2 $ENODE,l SEMODIF^ALL 

212. NALL SEALL 

213. /PNURMAT.l SLUMBER, 1 

214. /COM, 

215. /COM, FIND R THAT DIVIDES THE AREA FOR THE UPPER AND 

216. /COM, LOWER WINDING INTO EQUAL AREAS 

217. /COM, 

218. T1 - (R5*R5) ST2 - (R4*R4) ST3 - ((Tl+T2y2) 

219. R - SQRT(T3) SDR45 - (R5 - R4yNR45 SD452 - DR45/2 

220. /COM, ***** 

221. /COM, THT- THETA WIDTH OF A TOOTH (NOT INCLUDING SHOE) 

222. /COM, THC-THETA WIDTH OF CONDUCTOR 

223. /COM, THT - 4.331 STHC - 6.791 

224. THT-4STHC-6 

225. /COM, 

226. /COM, INSERT THE A+ MATERIAL IN THE BOTTOM OF SLOTS 1 & 2 

227. /COM, THE C-MATERIAL IN THE BOTTOM OF SLOTS 3 &4 

228. /COM, AND THE B+ MATERIAL IN THE BOTTOM OF SLOTS 5 & 6 

229. /COM£LOT#l 

230. MAT*5 

231. Tl-THT/2 ST2 - T1 + THC + EPS ST1 - T1 - EPS 

232. NSEL*X(R - EPS),(R5 + EPS) SNRSEL, Y,T1,T2 

233. ENODE,l SEMODIF r ALL 

234. NALL SEALL 

235. /COM^LOT #2 

236. T1« T2 + THT ■ 2*EPS ST2 - T1 + THC + 2*EPS 

237. NSEUX,(R - EPS),(R5+EPS) SNRSEL, Y.TLT2 

238. SENODE.1 SEMODIFrALL 

239. NALL SEALL 

240. /COM, SLOT #3 

241. MAT.10 

242. T1 - T2 + THT - 2*EPS ST2 - T1 + THC + 2*EPS 

243. NSEUX,(R - EPS),(R5+EPS) SNRSEL,Y,T1,T2 

244. ENODE,l SEMODIFrALL 

245. NALL SEALL 

246. /COM^LOT #4 

247. T1 - T2 + THT - 2*EPS ST2 - T1 + THC + 2*EPS 

248. NSEUX,(R - EPS),(R5+EPS) SNRSEL,Y,T1,T2 

249. ENODE,l SEMODIFrALL 

250. NALL SEALL 

251. /COM, SLOT #5 

252. MAT,7 

253. T1 - T2 + THT - 2*EPS JT2 - T1 + THC + 2*EPS 

254. NSELX/R - EPS),(R5+EPS) SNRSEL»Y,T1,T2 

255. ENODE.1 SEMODIFrALL 

256. NALL SEALL 

257. /COHSLOT46 

258. T1-T2 + THT-2*EPS ST2-T1 + THC + 2*EPS 

259. NSEUX^R - EPSX(R5+EPS) SNRSEL,Y,T1,T2 

260. ENODE,l SEMODIFRALL 

261. NALL SEALL 

262. /COM, 

263. /COM, NOW INSERT THE A+ MATERIAL IN THE TOP OF SLOTS 1 & 2 

264. /COM, THE C- MATERIAL IN THE TOP OF SLOTS 3 Sc 4 

265. /COM, AND THE B+MATERIAL IN THE TOP OF SLOTS 5 &6 

266. /COM, NEW MODIFICATIONS ON 1/26/93 INDICATE FULL PITCH WINDINGS 

267. /COM.SLOT #2 

268. MAT,5 

269. T1 - THT* 1.5+THC-EPS ST2-T1 + THC + 2*EPS 

270. NSEL^(R4-EPS)^R+D452) SNRSELY.T1.T2 

271. ENODE.1 SEMODIFrALL 

272. NALL SEALL 

273. /COM^LOT #3 

274. MAT, 10 

275. T1» T2 + THT - 2*EPS ST2 - T1 + THC + 2*EPS 

276. NSELX(R4-EPS),(R+D452) SNRSEL,Y,T1,T2 

277. ENODE,l SEMODIFrALL 

278. NALL SEALL 

279. /COM^LOT #4 
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280. MAT,10 

281. T1 “ T2 + THr - 2*EPS $T2 - T1 + THC + 2*EPS 

282. NSELX^R4-EPS),(R+D452) SNRSEL,Y,T1,T2 

283. ENODE.1 SEMODIF,ALL 

284. NALL SEALL 
283. /COM^SLOT #3 

286. MAT,7 

287. T1 - T2 + THT - 2*EPS $T2 - T1 + THC + 2*EPS 

288. NSELX(R4-EPSX(R+D452)SNRSEL > Y f Tl f T2 

289. ENODE.2 SEMODIF^ALL 

290. NALL SEALL 

291. /CORSLOT#6 

292. MAT,7 

293. T1 - T2 + THT - 2*EPS ST2 - T1 + THC + 2*EPS 

294. NSEL*X,(R4-EPS),(R+D452) SNRSEL,Y,T1,T2 
293. SENODE.1 SEMODIF^ALL 

296. /COM£LOT #1 

297. NALL SEALL 

298. MAT.5 

299. T1 - THT/2 - EPS ST2 - T1 + THC + 2*EPS 

300. NSEL^(R4-EPSX(R+D452)SNRSEL,Y,T1,T2 

301. SENODE,l SEMODIF,ALL 

302. NALL SEALL 

303. /PNUM,MAT,1 S/COLOR, 1,CYAN S/COLOR^BLUE S/COLOIU,R£D 

304. /COLOR,4,MAGENTA S/COLOR,5,MRED S/C0L0R,63MAGENT 

305. /COLOR,7,CBLUE S/COLOR,8,GCYAN S/COLOR,9,GREEN S/COLOR, 10, GREEN 

306. /NUMBER, 1 SNALL SEALL 

307. /COM, 

308. /COM, ENTER BOUNDARY CONDITIONS 

309. /COM, 

310. NSEI^XJU) 

311. NT,ALL»MAG,0 

312. NALL SEALL 

313. NSEUX36 

314. NT,ALI*MAG,0 

315. NALL SEALL 

316. CE,lAl>(AO,l,l^l>lAG,l 

317. RP30,1„U1 

318. /COM, 

319. /COM, ENTER THE CURRENTS 

320. /COM, 

321. JAPL-8.E6 

322. JBPL- -4.0E6 

323. JCPL - -4.0E6 

324. JAM • -JAPL SJBM--JBPL SJCM--JCPL 
323. KTEMP,-1 

326. NALL SEALL 

327. ESEL*MAT4 STE^ALL^JAPL 

328. NALL SEALL 

329. ESEUMAT.6 STE^LL^JBPL S/COM^HOULD BE 

330. NALL SEALL 

331. ESEL*MAT,7 STE^LUJCPL 

332. NALL SEALL 

333. ESEUMAT,8 STE*ALL*JAM 

334. NALL SEALL 

335. ESEUMAT^ STE^LUJBM 

336. NALL SEALL 

337. ESEUMAT,10 STE^LI^JCM 

338. NALL SEALL 

339. KBC,0 

340. TTER,10 

341. LWRTTE 

342. ITElL-20^2 

343. CNVR,0.0001 SKBC.1 

344. LWRITE SAFWRI 

345. FINI 

346. /EXEC 

347. /INPUTS 

348. FINI 

349. /EOF 
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APPENDIX G 


Experimental Thermal Investigation 


To start to become familiar with techniques and problems to be encountered in future 
laboratory work a commercially available motor was obtained and subsequently opened i.e. 
the end bells removed. From this first phase it was observed that there were eight poles and 
twenty seven slots. Thus it is immediately clear that this is a fractional slot motor i.e. there 
is no integer by which the product of the number of phases and the number of poles will 
equal 27. It was interesting to see that the manufacturer had used extended end windings on 
one end of the motor and carefully formed them so that they fit tightly into the aluminum 
end bell (without however any potting). Based on our later work this is, no doubt, to 
facilitate heat removal. Next the end windings on one end were cut off more or less flush 
with the first laminations in the stator stack. A 600 dpi color scan was made of the end of 
the lamination stack and is enclosed as figure on p. 99 below. The numbers beside some of 
the slots denote the numbers of conductors in the various slots. From their varying numbers 
it is clear, as expected, that this fractional slot motor does not have identical coils. Thus its 
electrical balance is assured from careful design of the unbalanced windings. The diameter 
of the wire was measured several times and it approximates 25 AWG but not exactly. 
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Perhaps the wire is si metric size (this is s pure speculation). Three Hall effect devices were 
found to be cemented to the inside of the assembly where the magnetic flux from the rotor’s 
permanent magne ts can sweep over them. The angular dis t anc e between them corresponded 
to the standard required for electronic commutation purposes. Additional things to be 
noticed the fact that there is no apparent effort made to encapsulate the wires in the 

slots. This is a little surprising because it is common commercial practice to vacuum 
impregnate the windings with insulating varnish after they are placed in the stator. In this 
inctanry* where heat removal was pretty clearly a design feature the voids in the slots is a 
little p uzzlin g (or perhaps it is felt that the large end turns and so forth are enough of a heat 
removal design feature to make the motor competitive). Analysis of the fuller slots showed 
that the fill factor achieved in this motor does not exceed 50% (a common number it 
appears). The stator laminati ons measured to be 25 mils thick. This is a little surprising 
because 14 mil material is generally available (10 mil material appears to be available in a 
few cases). However in the light of our core loss analysis perhaps core loss is just not that 
big a problem and possibly the thicker material is more easily fabricated than the thinner. 

After this examination it was decided to perform a simple thermal test on the motor. The 
motor was fixed to the bottom of a plastic jar with RTV. Thermocouples were cemented 
along the length of the inside of the stator and to the copper end windings (black and white 
photographs). The first thermocouple was placed 1/4 inch up from the bottom of the stack, 
the second 3/8, third 3/4, fourth 1 1/8, fifth 1 1/2, and the last two on the copper at the top of 
the motor. With a known ambient temperature as an initial condition an ice and water 
combination was poured into the plastic jar Le. surrounding the outside of the motor and the 
output of the thermocouples recorded on a floppy disk. The results of this exercise are 
presented in the mul tiline multicolored figure on p. 101 showing temperature at the various 
stations versus time. 
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While the exponential appearing changes in the various temperatures confirms the nature of 
the response and agrees with the lumped thermal circuit approach to modeling one thing 
stands out that is somewhat counterintuitive (at least to this author without the benefit of a 
ther mal simulation of this configuration). That is die temperature of the copper followed the 
temperature of the iron without any appreciable time lag. While there is at present no 
explanation of this observed fact speculation must center on two things. The first 
speculation is that there is better than anticipated iron to copper wire thermal coupling 
(despite the visible voids in the slots). The second speculation is that the relatively small 
thermal capacitance of the copper (compared to the relatively massive iron presence) makes 
it possible for its temperature to track the iron temperature closely in time. 

Several possibilities for instrumenting the fabricated motors arise from this experiment. To 
perform meaningful thermal characterization of the new motors temperature measurements 
must be made in several contexts. Of course thermal measurement wires should be wound 
with the main windings (this need only require a single strand). During laboratory tests 
essentially the same experiment described above should be run with refinements gained 
from this work. One experiment might involve thermocouples cemented to the iron and the 
copper as before. However instead of not only cooling the outer part of the stator and seeing 
how heat leaves the motor in an attempt to characterize the thermal circuit but also the 
copper winding would be energized to determine its temperature profile with time. Another 
test to be done is to spin the unenergized motor with another motor and measure its 
temperature rise and, if possible, the power transmitted into the motor. This test would 
establish the core loss caused by spinning the permanent magnets relative to the stator and 
confirm hopefully that it is a small contributor. 
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Appendix H 


Motor Magnet Restraint 

During this Delivery Order Professor Campbell taught a senior level design course in 
engineering. The interest of one of his students, Mr. Daniel DeHaye, was piqued by the 
problem of shrink fitting thin, metallic, non megnetic tubing over the rotor magnets to 
prevent them coming loose from the rotor especially during high speed motor operation 
(where failures have been reported occuring in other designs). Thus Mr. DeHaye pursued an 
investigation of this problem (albeit the 4 pole design) for his class project 

His final report is presented in this appendix. Only the color figures have been added to his 
original work to aid in interpretation. The ANSYS system used for this analysis is a 
dimension limi ted one resident on a Hewlett Packard machine at UAH. It twill be seen that 
he has solved the problems of not only the stresses and strains but also of including the 
effects of rotation and the temperature to which it is required to heat the sleeve prior to 
placing it over the rotor. His results show this is a very feasible approach to magnet 
retention. 

Mr. DeHaye's figure of 15-20% cost of Shuttle refurbishment due to hydraulics is too high. 
However it is substantial and a NASA effort to cut or eliminate it is actively underway. 
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Motor Magnet Restraint 


ME 465 
Final Report 
3-12-93 


David DeHaye 



I. Problem description 


This project is part of an effort by the Research Institute (RI) and the University 
of Alabama in Huntsville (UAH) to design a 3 inch permanent magnet motor. The 
objectives of my project were to define a methodology for the design of magnet restraints to 
be used on permanent magnet motors, and to use this methodology to design magnet 
restraints for the RI/UAH motor. The Rl/UAH motor is designed to be used on an 
electromechanical thrust vector control actuator. The electromechanical system is intended 
to replace hydraulic systems that can be expensive to maintain. For example, each time a 
space shuttle is refurbished, 15-20% of the cost comes from the hydraulic systems. 
Therefore, my project will contribute to a design effort that could save a large amount of 
money. The magnet restraint is being designed to overcome problems with motors that use 
magnets that are bonded to the rotor using adhesives. The problem with using a motor that 
has glued on magnets is that at high rpm the centrifugal forces involved can cause the 
magnets to break free from the rotor at the bond line. The restraint being designed is a non¬ 
magnetic sleeve that will be shrink fitted over the magnets to hold them onto the rotor. By 
using finite element modeling (FEM), a shrink fit can be designed so that the magnets are 
always in compression at their bases. This will ensure that the magnets will not break free 
from the rotor. 

n. Design options considered 

When I started to work on this project, the method of restraint had already been 
decided. A non-magnetic sleeve will be shrunk over the magnets to hold them in place. 
However, the method of performing the design was not specified. There are at least three 
methods that could have been used to perform the design. It could have been done by hand, 
but this method might give inaccurate results and it would require that the entire process be 
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repeated for different motor configurations, such as a larger motor or a different number of 
magnets. The second method is to write a fill in the blanks computer program. This would 
be an automation of the process done by hand. The accuracy of the results would be no 
better, but once the p rogram is written, the analysis could be repeated quickly for different 
motor configurations. The third method is to use an FEM computer program such as 
ANSYS. This method should provide accurate results and after the model is first defined, 
only a few lines of the input file will need to be changed to account for different variables. I 
chose to use ANSYS to do the design of the motor magnet restraint 

TTT Reasons for selection of design method 

I chose to do the design of the motor magnet restraint using ANSYS. There are a 
number of reasons why I chose this design method. I am familiar with the program, and the 
results are probably more accurate than the results I would get using other methods. Also, 
different motor configurations can be analyzed by changing a few lines in the input file. 
Another reason is that the shrink fit is easy to simulate using ANSYS. This is done by using 
an interface element that can specify an interference fit between two adjacent elements. 
These reasons lead me to believe that the use of ANSYS is the best method for doing the 
design. 

TV. Design description 
A. Design overview 

Figure I shows what the cross section of a typical permanent magnet motor looks 
like. The magnets are glued to the rotor and can break loose at high rpm. This problem is 
going to be fixed by designing retaining sleeves that will be shrunk over the magnets to 
hold them onto the rotor. An element plot of the restraint and motor can be seen in figure 
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2. The bands of thin elements at the outside edge of the circle represent the retaining 
sleeve. The rest of the elements represent the magnets, the rotor, and the motor shaft 

The design of the retaining sleeve for the 3 inch Rl/UAH motor was done using 
ANSYS. Titanium was used for the sleeve because it is light in weight relatively stiff, 
non-magnctic, and has a high yield strength. I found that by using an interference fit of 
one part in a thousand, and a sleeve thickness of 1 mm, the magnets would stay in 
compression at their bases at 20,000 rpm. This setup gives a factor of safety of 6.5. This 
amoun t of interference will require that the sleeve be heated to about 110°C (200°F) 
above the temperature of the motor to do the shrink fit This should be acceptable if the 
sleeve is always at the same temperature as the motor and magnets. If the temperature of 
the sleeve should raise above the temperature of the motor, the shrink fit could loosen. If 
this is a potential problem, the shrink fit could be made with a greater interference so 
that at operating temperatures it keeps the magnets in compression at their bases. I ran 
the problem using a 1 mm thick ti tanium sleeve and the interference that comes from a 
210°C (380°F) shrink fit The factor of safety for this configuration was 3.4. For both of 
these configurations, the magnets remain in compression at their bases, the factor of 
safety is good, the sleeve is adequately thin, and the shrink fit temperatures are 
reasonable. Therefore, a titanium sleeve of 1 mm in thickness should be adequate to 
retain the magnets on the 3 inch RI/UAH motor. Also, the use of ANSYS as the method 
for doing the design seemed to be successful. 

B. Task descriptions 

There were two tasks required to finish this project The first was to write an input 
file for ANSYS to do the design of the restraint The second was to run the program for 
various motor configurations and analyze the data that was generated by the program. 
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The input file for ANSYS and a description of the variables that can be modified is 
shown in appendix A. In writing the input file I encountered mostly minor problems. 
The biggest problem that I had was in coming up with a way to s imulate the interference 
fit between the sleeve and the motor. I discussed this problem with Dr. Campbell and he 
called the technical support people at ANSYS. They suggested that an interface element 
be used. This type of element is used between surfaces that can make and break contact 
with each other. The attractive feature of this element is that it allows for an interference 
to be specified between two adjacent elements. By placing these elements between the 
inside of the retaining sleeve and the outside of the magnets, a shrink fit can be 
simulated. To see if this type of element could be used to simulate a shrink fit, I tested it 
on a simple problem. The test problem was a press fit of a ring onto a shaft The results 
obtained using ANSYS and the interface element were within five percent of the 
textbook answer. This convinced me that interface elements could be used to simulate 
the shrink fit of the restraint over the motor and magnets. 

I ran the problem for several motor configurations and analyzed the output. The output is 
in the form of stress plots and stress tables. For each configuration the first step in 
analyzing the restraint system was to make sure that the magnets remained in 
compression at their bases. I did this by looking at the stresses for the nodes on the bases 
of the magnets. If the stresses in the radial direction are compressive for all of these 
nodes, then the magnets must be in compression at their bases. A listing of the stresses 
for the nodes at the base of the magnets in the 3 inch RI/UAH motor is given in 
appendix B. This listing is for a titanium sleeve that is 1 mm thick and has an 
interference of one part in a thousand. SX is the stress in the radial direction and the 
units are in pascals. Once it has been determined that the magnets are in compression at 
their bases, it is necessary to determine whether the sleeve can tolerate the stresses 
involved. For each configuration the maximum stress intensity is found by making a plot 
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of the stress intensity. A stress intensity plot for the 3 inch Rj/UAH motor is shown in 
figure 3. The maximum and minimum stress intensities are marked by an MX and MN 
on the plot To the side of the plot, the values for MX and MN are given in pascals. The 
factor of safety is determined by dividing the yield strength of the sleeve material by 

MX. In this case, titanium has a yield strength of 830 MPa. Therefore, the safety factor = 

830.0MPa , . .... . .. 

-- 6.5, which is acceptable. 

127.8MPa 

This design process was done for several other motor configurations. Some 
interesting results were obtained by doing this. For a motor that is roughly twice as big 
in diameter as the RI/UAH motor, the factor of safety using a I mm thick titanium sleeve 
was found to be 1.75. The interference fit had to be increased to about one part in 250 to 
keep the magnets in compression at their bases. This made the maximum stress intensity 
go up and the factor of safety go down. This also means that the sleeve would have to be 
about 420°C (750°F) above the temperature of the motor to do the shrink fit. These 
results indicate that the use of a shrink fitted retaining sleeve might be limited to smaller 
motors. I had one interesting problem in doing the analysis of the RJ/UAH motor. The 
man ufacturers data for the magnets gave a modulus of elasticity that was about one 
thousand times less than that of steel. I thought that this was too low and the company 
was called to find out whether a mistake had been made. They said that this was the 
correct value, so I ran the problem using this modulus of elasticity. A stress intensity 
plot of the motor and restraint showed that it had deformed into an almost square shape. 
This plot is shown in figure 4. Another call was made to the company and this time they 
said that the value given was a thousand times too low. This problem took about a week 
to resolve and it is a good example of why you should always check your data and facts. 
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C Design integration 

There were several constraints and requirements that had to be satisfied for the 
successful completion of this project, these are listed below. 

1) The material for the sleeve has to be non-magnetic. 

2) The sleeve must be as thin as possible so that the air gap of the motor is not 
affected. 

3) The magnets should remain in compression at their bases. 

4) The temperature required to perform the shrink fit should be low enough that the 
magnets will not be adversely affected. 

5) The method used for the design should be ada p ta b le to other motor 
configurations. 

6) The Rl/UAH motor should be able to withstand speeds up to 20,000 rpm. 

D. Design safety 

If the retaining sleeve were to fail in service it would cause the electromechanical 
thrust vector control actuator to stop working. This could result in a scrubbed flight or 
possibly even the crash of the vehicle that it was used in. In either case, the monetary 
cost would be very high. If the vehicle crashed, this could cause injuries or death. 
However, the design will probably go through a great deal of testing before it will ever 
be used on an actual launch vehicle. Therefore, if the retaining sleeve does fail, the costs 
will be limited to the cost of a prototype and the cost of testing. 
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V. Summary 


I believe that the two main objectives of this project have been met. The first 
objective was to define a methodology for the design of magnet restraints to be used on 
permanent magnet motors. I did this by using ANSYS to do the design. It gives accurate 
results and it is easy to handle different motor configurations by changing the input file of 
the program. The second objective was to use this methodology to do a design of a magnet 
restraint fen- the 3 inch RI/UAH motor. I used ANSYS to design a restraint for this motor 
and determined that a I mm thick titanium sleeve will work. Therefore, I think that this 
project was a success. If I were to do it over again, I would not do it differently, but I could 
have used more time. With extra time I could have investigated the possibility of using a 
thinner sleeve. 
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Motor Cross Section 
Showing Finite Elements 
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Rotor Stress Plot 
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Enlarged View of Stresses 
In Titanium Sleeve 
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MOTOR MAGNET RESTRAINT ! 
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Program listing 
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nafgr/'g^ p roterf/es. 

E. ^ c7 cfasTic/T^ 

OEMS- c/<r/7 <; 7 y 


hi \jXy - Pc> <:• Sf<< rpTr'c 
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***7V/■<■'*/ ij 


>~toTer jUaft 
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r+u~f ri£ / % 



tUe r**cT c r ;li*ft. 


1~he mcTot T 


■Fro/*? ile ce^l/dT i^tc .^>otcr 

f~o -H*&. ciJTs/d^. of The. /~c~Tc 




rfrdiuf the carver He 

To fite CxJTs'J* £>F *6* r*r*r*?£T; 






Sane a£ ahe?is£ p!\*.i a J*$Teir 

r&-dsit£ ~f~t'c* rr > -the cerVcf cT rite stycTt 
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Appendix B. 

Nodal stresses at magnet bases 


140 


P0ST1 NODAL STRESS LISTING 


LOAO STEP 1 ITERATION- 10 SECTION- 1 

TINE- .00000E+00 LOAD CACC- 1 

THE FOLLOWING X,Y,Z STRESSES ARE IN COORDINATE SYSTEM 1 


NODE SX SY SZ SXY 

SIG2 SIG3 SI SIGE 

235 - . 6214E* 07 -.5086E+07 -.318DC-*07 .82621:-<06 

L T >-07 - . 3555E+07 -.7832E»07 .4331E+07 .4604E+07 

240 - . 621 0E ♦ 07 -.5034E*07 -.31B6E«-07 -.3232E+06 

L<07 - . 3553E* 07 -,7878E-*07 .4829E<07 .4602C<07 

245 -.4694E*07 -.2517E*07 -.2016E+07 -,24?2E*06 

L<07 -.2440E + 07 -,4771E-<07 .2755T" 07 .2574C"07 

250 -.3256E* 07 -.2202E+07 -.1525E*07 ,1995E*D6 

'••'07 - . 2129E ♦ 07 -.3329E-*07 ,1C04E-<07 .1592C-<07 

255 - . 269 IE* 07 -.2036E*07 -.1321E*07 .1162E+06 

f>07 - . 1980E* 07 - . 2747C-* 07 .1426E+07 .1238C<07 

260 -.2553E+07 -.1936E*07 -.1267E*07 4093. 

f-:^07 -.1957E+07 -.2581E+07 .1314T107 .1146f>07 

265 - . 2685E + 07 -.2050E*07 -.1323E>07 -.1135^1-06 

07 - . 1991E+07 - . 2744C + 07 .1421C<07 .12330^7 

270 -.3262E*07 -.2220E*07 -.1532E»07 -.2063E*06 

L • 07 - . 2142E + 07 -.3339E«07 .1808C<07 . 1594C-< 07 

275 -.4725E+07 -.2534E+07 -.2029E+07 .2356E+06 

E<07 - . 2458C+ 07 -.4001E*07 .2771E-<07 .25GCCi07 

230 -.6273E+07 -.5098E+07 -.3209E+07 .8331E+06 

L<07 -.3586E+07 -.7930E+07 .4861E->07 .4629Et07 

MAX IMUMS 

NODE 280 280 280 280 

280 280 280 230 
UALUE -. 6278E* 07 -.5098E-<07 -.3209E+07 .P331E<06 

9E + 07 -.3536E + 07 -,7930E*07 .4361E+07' .4629E+07 


SY7 


SX7 


nr 

-.1349C- 

-05 

. 49421"- 

-n? 

-. 5 r> 

. 7253E- 

•1 7 

.1197E- 

1 9 

- . 50 

. 3408E- 

06 

. 7 76 3 F - 

ns 

- , ?n 

. 2663E- 

06 - 

- . 23-03E- 

ns 


-.1511E- 

06 

.17R3E- 

ns 

-.15 

.8751E- 

09 - 

-.2743E- 

0 6 

- . 1? 

.ooooe+ 

00 

. O000E + 

no 

- . 15 

.3106E- 

06 

.2599E- 

ns 

- . 1 S 

.4540E- 

06 - 

■.5193E- 

05 

-,?n 

-.1423E- 

05 

.5243E- 

05 

- . 3 0 


190 190 2 

-.1472E-05 .5397E-05 -.3 


ORIGINAL PAGE !$ 
OF POOR QUALITY 
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P0CT1 NODAL STRESS LISTINC *»«** 


LOAD STEP 1 ITERATION- 10 SECTION- 1 

TIME- .OOOOOC+OO LOAD CASE- 1 


T!!E FOLLOUING X,Y,Z STRESSES ARE IN COORDINATE SYSTFM 1 


E 


tl 


NODE 

S 

170 
+ 07 
175 
+ 07 -. 
100 
♦07 
1C5 
+ 07 
190 
+ 07 
195 
-07 
200 

• 07 

205 

♦ 07 

210 
+ 07 
215 
+ 07 - . 
220 
+ 07 
225 
+ 07 
230 
+ 07 


SX SY S7 

IG2 SIG3 SI 

-.2553E+ 07 -.1976E+07 -.1265E+ 
1949E+07 -.25G1E+07 .1316E+07 

- . 2691E+ 07 -.2033C+07 -.1320C+ 

1978E+07 -.2746E+07 .1426E+07 

-.3257E+07 - . 22 01C+ 07 -.1525E+ 
2129E + 07 -.3329E + 07 .1804E + 07 

- , 4695E + 07 - . 25 18C + 07 -.2017E< 

2441E + 07 -.4772E + 07 .2756E+07 

-.6215E+ 07 -.50C9C+ 07 -.3189C+ 
3556E+07 -.7885E+07 .4333E+07 

-.6215E*07 -.5089E+ 07 -.3189C+ 
3557E+07 -.7885E+07 .4833E+07 

- . 4695E ♦ 07 - . 251 8C •< 07 -.2017E< 

2441E + 07 -.4772E + 07 .2756E+07 

-.3257E+ 07 -.2201C+ 07 -.1525C+ 
2129E + 07 -.3329E + 07 .1804E + 07 

-.2691E+ 07 -.2033E+ 07 -.1320E+ 
1978E+07 -.2746E + 07 .1426E+07 

-.2553E*07 -.1976E+ 07 -.1265Ci 
1948E + 07 -.25G1E + 07 .1316E+07 

-.2691E+07 -.2033E+ 07 -.1320C+ 
1977E+07 -.2746E+07 .1426E+07 

-.3256E + 07 -.2200C+ 07 -.1525E + 
2129E + 07 -.3323E + 07 .1804E + 07 

-.4694E+ 07 -.2517E+ 07 -.2016E+ 
2440E+07 -.4772E+07 .2755E+07 


SXY 

SIGE 

07 46.75 

.1147E+07 
07 -.1124C+06 
.1230E+07 
07 -.196CE+06 
.1592E+07 
0? . 2447E-! 06 

.2574E+07 
07 .C254T *06 

.4 6 06E + 07 
07 -.C256E+ 06 
.4606E♦07 
07 -.244CC+ 06 
. 2574E + fl7 
07 .1967F +06 

.1592E+C7 
07 . 1 1 ?4C-i 06 

.1230E+07 
07 59.76 

.1147E+07 
07 -.1123E+06 
.1238E+07 
07 -.1965C+ 06 
.1592E+07 
07 .2450C+06 

.2574E+07 


SYZ SX7 CIG 


- .1607C-10 

.4429E-06 


. 1 ? 

.48 08C-07 

.5863E-06 

- 

. Y 

.2957E-06 

. 2602C-05 

- 

. 1 f - 

. 4355C-06 

- . 4767F7- 05 

- 

. ?r 

-.1472E-05 

.5397r-05 

- 

1 t r 

.129 OC-16 

.4779C-16 

- 

. •* r 

.4367F-06 

.4777E-05 

- 

n f 

. 35 02C - 06 

- . 3003C-05 

- 

. 1 r 

-. 1609C-06 

.2060C-05 

- 

. 1 “ 

-. 2556E-10 

.551 IE-06 

- 

. Y 

. 3G58C-07 

. 4712C-06 

- 

. V 

. 2556E-06 

. 2253C-05 

- 

. <1* 

. 4721E-06 

-.5160C-05 

- 



OF POOR QUALITY 
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08 13:08.1991 MOUT.DAT Page 1 



:-'PlNT ALL v NODAL STRESSES PER NODE 

• < • 

.- #**** POSTl NODAL STRESS LISTING ***** 

LOAD STEP 1 ITERATION- 10 SECTION- 1 

1 I ME- ? iOOOOOC+OO LOAD CASE- 1 

T!IE following' x,v,z stresses are in coordinate systfm 1 


NODE SX SY SZ SXY 

SIG2 SIG3 SI SIGE 

105 -.61600*07 -.50800* 07 -.3171E+07 -.81600*06 

:■ i 07 -.3530E*07 -.>845E*07 .4303E*07 .4586E*07 

110 -.4669E* 07 -.25050*07 -.20060*07 -.24970*06 

: P.7 - . 2427E* 07 -,4748E*07 .2742E*07 .2561E*07 

115 -.32500+07 -.21080*07 -.1519E<07 .18990*06 

: :!j7 - . 2120E *-07 - . 33190*07 .1799E’07 .15800*07 

120 -.2691E* 07 -.2023E+07 -.13170*07 .1100C*06 

i_. : U7 - . 1970E*07 - . 2744E *07 . 1427E*07 .12390*07 

125 -.25550*07 -.19690*07 -.12630*07 649.6 

07 - . 1943E*07 -.2531E+07 . 131OE*07'' .1149E*07 

130 -.2692E* 07 -.20300*07 -.13190*07 -.11070*06 

L :u7 - . 1975E*07 -.27470+07 .1427E+07 . 12390*07 

135 - . 32580+07 -.22000+07 -.15250*07 -.19500*06 

- 07 -.2129E+07 -.33290+07 .1304E+07 .1592E+07 

140 - . 4696E* 07 - . 2520C+ 07 -.20170+07 .24640*06 

= 07 - .2442E + 07 -.4774E + 07 .2757E+07 .25750+07 

145 -.62190+ 07 -.5093E+07 -.3191E+07 .82720*06 

; m)7 - . 3559E *07 -.7890E+07 .4336E+07 .4603E+07 

150 - . 62150+ 07 -.5092E+07 -.3190C*07 - . 8249C* 06 

EhJ7 - . 3557E+07 -.7836E+07 .4834E+07 .46070+07 

155 -.4695C+ 07 -.25190*07 -.20170*07 -.24460*06 

;.: 07 - . 2441E *07 - . 4773E *07 .27560+07 .2574E+07 

160 -.3257E+ 07 -.22010+0? -.1525E+07 .19690+06 

i; :■ 07 -.2129E*07 -.3329E+07 .1804E+07 .15920+07 

165 -.26910+07 -.2033E+07 -.13200*07 .11250*06 

l:-07 - . 1973E *07 -.2746E+07 .1426E*07 .12330*07 


SYZ 

SXZ SIR 

.14160-05 

.51050-05 -.30 

.42650-06 

.45500-05 -.20 

.36590-06 

-.33350-05 -.15 

-.16530-06 

.20640-05 -.13 

.00000*00 

.00000*00 -.12 

-.23600-07 

-.29350-06 -.13 

.25370-06 

.2254T- 05 -.15 

.32C6C-06 

-.35720-05 -.20 

-.11640-05 

.42600-05 -.30 

.11440-18 

.41940-18 -.30 

.39940-06 

.43740-05 -.20 

.37940-06 

-.33370-05 -.15 

-,14640-06 

.17040-05 -.13 


ORIGINAL PAGE IS 

OF POOR QUALITY 
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APPENDIX I 


Inertia Considerations 


Mention has already been made in the main body of the report about inertia concerns. For 
ready refemece are listed immediately below the results of optimizing the motor, spur gear, 
roller screw configuration for maximum acceleration. 


U1 -> 


4 J« Sqrt[Jr») 

Sqrt (-+ 4 SqrtfJp) Sqrt(Jrs) + 8 Jrs] PI 

Sqrt[Jp) 

Sqrt(J1) 


1/4 

Jr« 

n ->-j, a -> 

1/4 

JP 


The various Js in the formulas represent the various inertis and reflected masses in the 
system. The spur gear ratio is denoted by n and the roller screw ratio by 1. See the original 
work for all the definitions. These results make it plain that die interrelationships are subtle. 

The following material exerpted from the previous report is included to lend some heuristic 
insight into the foregoing relationships. 
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It has been stated elsewhere in this report that to obtain maximum possible performance 
from a motor, gear train and load combination it is necessary to match carefully the motor to 
the load by careful selection of a gear ratio between the two and then repeat the matching 
process for a number of differing motors to ascertain the best possible load acceleration 
performance. The results may be counter intuitive to a casual observer because in many 
minds it is natural to think in terms of the motor and its speed and acceleration rather than 
the speed and acceleration of the load. In the present design environment it is possible to 
select motors of the same rated horsepower with different available torque capacities 
because of different speeds at which the motors are rated (and they therefore have differing 
attendant polar mass moments of inertia) as well as being able to select independently the 
gear ratio(s). Thus in the present case only the load is a fixed quantity - and it is the load in 
whose motion one is intrerested. The results of this procedure are presented in this report 
for the multipass gear train case but a simpler example will perhaps help to cement the 
principles involved. 

Consider the case of a motor coupled to a load by a simple one pass gear train. Now 
examine the accelerations of the motor and the load. The system inertia referred to the 
motor shaft will be given by 


Jtm=Jm+Jl/n~2 

while the system inertia referred to the load is given by 

Jtl=n~2 Jm+Jl 
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Where n is a step down gear ratio from the motor to the load and is greater than one in 
magnitude, Jt denotes a total inertia either at the load (1) or the motor (m), JM is the motor 
rotor inertia and J1 is the load inertia. 

The expression for the motor acceleration is given by 

accel(m)=T/Jtm=T/(Jm+Jl/n“2) 


and the load acceleration is 


accel(l)=Tn/(Jl+Jmn*2) 

where T is the torque available from the motor. 

Inspection of these expressions show immediately that to maximize the acceleration of a 
given motor n should be made as large as possible; in the limit the motor becomes 
uncoupled from the load (Jl/n~2 approaches 0). Finding the maximum acceleration of the 
load is a more tedious matter but if the derivative of the load acceleration is taken with 
respect to n and the results set to zero the results become 

n=SquareRoot[JLJm] 

accel(l)max= T/(2 SquareRoot[Jl Jm]). 

Thus even in the simplest case finding the correct match of inertias through gear train ratio 
selection is more involved than maximizing the motor acceleration. And in this case for a 
given load different motors with different maximum torques and different inertias will 


146 



produce different maximum load acceleration (which is desired for maximum bandwidth or 
speed of response). In this maximized load acceleration case the acceleration of the motor 
will be (by simple substitution) expressed by 

accel(m)=T/(2 Jm) 

which is certainly lower than could be achieved by a higher gear ratio (as n approaches 
infinity accel(m) approaches T/Jm). 

It is hoped that the simple example discussed above will help the intuitive reasoning process 
concerning motor and gear ratio selection for a given load. Hopefully it lends insight into 
why the multipass gear case is somewhat involved to optimize with respect to load 
acceleration and motor-gear train selection, as has been done in this report. 
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